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A  COMPUTER  SOLUTION  FOR  THE  DYNAMIC  LOAD, 
LUBRICANT  FILM  THICKNESS  AND  SURFACE 


TEMPERATURES  IN  SPIRAL  BEVEL  GEARS 


SUMMARY 


A  complete  analysis  of  spiral  bevel  gear  sets  Is  presented. 
The  gear  profile  Is  described  by  the  movements  of  the  cutting  tools. 
The  contact  patterns  of  the  rigid  body  gears  are  investigated.  The 
tooth  dynamic  force  is  studied  by  combining  the  effects  of  variable 
teeth  meshing  stiffness,  speed,  damping  and  bearing  stiffness.  The 
lubrication  performance  is  also  accomplished  by  including  the  effects 
of  the  lubricant  viscosity,  ambient  temperature  and  gear  speed.  A 
set  of  numerical  results  is  also  presented.  s' 
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CHAPTER  I 


INTRODUCTION 

Spiral  bevel  gears,  found  in  many  machine  tools,  automobile 
differential  gears, and  aircraft  gears,  are  important  elements  for  power 
transmissions.  However,  the  basic  mechanisms  which  govern  the  major 
failure  modes  of  spiral  bevel  gears  are  not  fully  understood.  Unlike 
spur  or  helical  gears,  the  complicated  geometry  of  spiral  bevel  gears 
makes  the  problem  considerably  more  difficult.  In  military  applications 
such  as  gear  transmissions  used  in  recent  V/STOL  aircraft,  gears 
are  often  designed  under  conditions  very  close  to  the  failure  limits 
to  gain  a  maximum  power  density.  A  thorough  understanding  of  spiral 
bevel  gears  under  critical  operation  is  urgently  needed  to  prevent 
premature  failure. 

Basically,  there  are  two  classifications  of  gear 
failure,  one  of  which  is  structural  failure,  which  includes  flexure 
fatigue,  tooth  breakage  and  case  crushing,  and  the  other  .is  the 
lubrication  related  failure  which  includes  wear,  surface  pitting  and 
scuffing. 

Structural  failures  are  usually  related  to  the  material  properties 
and  improper  geometrical  design  or  unexpected  overloading.  While  the  relia¬ 
bility  of  spira.1.  bevel  gears  can  be  improved  by  better  material 
selection  and  better  geometrical  design,  it  can  also  be  enhanced  by 
a  more  accurate  understanding  of  the  dynamic  loading  between  a  pair  of 
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teeth  as  it  travels  through  the  contact. 

Since  the  basic  lubrication  process  in  bevel  gears  is  not  clearly 
understood,  the  lubrication  related  failures  .particularly  in  surface 
pitting  and  scuffing, are  much  more  difficult  to  predict.  Current  methods 
used  in  industries  to  predict  lubrication  related  failures  are  not 
satisfactory  because  they  are  mainly  empirical.  It  has  been  shown  in 
failure  tests  of  gears  and  rollers  that  the  surface  pitting  and 
scuffing  are  affected  critically  by  lubricant  film  thickness  and  surface 
temperature.  A  tool  that  can  predict  film  thickness  and  surface  tempera¬ 
ture  in  the  gear  teeth  contact  is  welcome  in  design  against  pitting  and 
scuffing. 

This  paper  presents  an  analysis  to  predict  the  dynamic  loads 
between  gear  teeth  contact  which  has  a  contact  ratio  greater  than  unity 
and  a  variable  stiffness  along  the  path  of  contact.  In  addition,  an 
analysis  is  presented  to  predict  the  minimum  film  thickness  and 
the  surface  temperature  under  given  dynamic  loads.  It  is  based  on  the 
most  recent  theories  on  film  thickness  and  traction  in  elastohydrodynamic 
contacts. 


CHAPTER  II 


TOOTH  PROFILE  AND  CONTACT  PATTERN 
OF  SPIRAL  BEVEL  GEARS 

2.1 _ Introduction 

Recent  developments  in  the  finite  element  method  (FEM)  and 
elastohydrodynamic  lubrication  (EHD)  have  stimulated  a  growing  interest 
in  the  basic  research  of  the  dynamics  and  lubrication  process  in  spur 
and  helical  gears  (1-6).  This  interest  is  now  gradually  being  extended 
to  spiral  bevel  gears  for  which  there  seems  to  be  little  known  about 
the  tooth  load  sharing  and  lubrication  characteristics,  especially  at 
high  speeds  and  heavy  loads. 

Unlike  involute  spur  and  helical  gears  whose  surface  geometry 
is  relatively  simple  and  readily  derivable,  the  surface  geometry  of 
spiral  bevel  gears  is  extremely  complex  owing  to  its  three  dimensional 
characteristics.  For  studying  the  performance  and  failure  of  spiral 
bevel  gears,  it  is  essential  first  to  find  an  efficient  method  to 
describe  the  surface  geometry  of  a  pair  of  mating  gears  in  terms  of  the 
tooth  form  and  the  tooth  contact  patterns.  Although  this  method  is 
currently  available  at  Gleason  Gear  Works  (7-11),  it  is  not  fully  suffi¬ 
cient  for  studying  the  dynamic  load  and  lubrication  analysis.  It  is  the 
intent  of  this  research  to  extend  the  existing  surface  geometry  analysis, 
known  as  the  Tooth  Contact  Analysis  (TCA),  to  generate  additional 
quantities  needed  for  the  lubrication  and  dynamic  analysis. 


The  formulae  for  obtaining  the  coordinates  of  tooth  surfaces 
based  on  the  Gleason  generating  system  a  re  derived  in  Section  2.2.  The 
position  of  each  point  on  the  surface  is  described  in  cylindrical 
coordinates  in  terms  of  its  position  along  the  rotational  axis  of  the 
pinion  or  gear,  its  radial  distance  from  the  axis,  and  its  angular  posi¬ 
tion  about  the  axis  from  an  arbitrarily  chosen  reference  plane.  The 
equations  are  solved  by  an  iterative  method. 

The  contact  points  of  a  pair  of  mating  gears  and  pinions  for  a 
running  position  are  determined  in  Section  2.3.  The  sliding  velocity 
and  the  direction  of  the  unit  normal  vector  at  each  contact  point  can 
be  obtained.  In  addition,  the  principal  radii  of  the  mating  surfaces  which 
characterize  the  geometry  at  the  contact  point  can  be  computed. 

Some  numerical  results  of  a  typical  case  are  shown  in  Section 
2.4.  The  effect  of  different  running  positions  on  contact  patterns 
are  studied.  Finite  element  grids  of  a  typical  spiral  bevel  tooth 
profile  are  presented. 
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2.2 _ Determination  of  Tooth  Surfaces  for 

Generated  Spiral  Bevel  Gears 

In  Tooth  Contact  Analysis  (TCA)  the  gear  and  pinion  tooth 
surfaces  are  mathematically  described  by  the  machine  settings  and  the 
cutter  specifications  used  in  making  the  gears.  Before  proceeding  to 
describe  TCA,  it  is  first  desirable  to  review  some  basic  concepts  of 
bevel  gearing.  The  principal  reference  planes  of  a  spiral  bevel  gearing 
system  are  shown  in  Fig.  2-1.  The  axial  plane  of  a  single  gear  is  the 
plane  containing  the  gear  axis  and  a  given  point  on  the  tooth  surface. 
The  pitch  plane  is  the  plane  perpendicular  to  the  axial  plane  and  tan¬ 
gent  to  the  pitch  circle.  The  transverse  plane  is  perpendicular  to  the 
axial  plane  and  the  pitch  plane.  Figure  2-2  and  Figure  2-3  show, 
respectively,  the  sectional  views  in  the  axial  and  transverse  plane  of 
a  bevel  and  pinion  system,  and  illustrate  the  most  common  terms  used  in 
bevel  gears. 

In  order  to  understand  the  formulae  derived  in  this  chapter,  a 
brief  explanation  of  the  cutting  process  is  in  order.  Figure  2-4 
illustrates  the  Gleason  cutting  process  —  1)  the  wheel  cutter  used  in 
cutting  spiral  gears,  2)  the  cradle  which  positions  the  cutter  in 
reference  to  the  work  and  carries  it  through  its  generating  motion, and 
3)  the  blank  into  which  the  teeth  are  cut.  The  dresser,  which  is  used 
to  restore  the  cutting  action  of  the  wheel  by  fracturing  and  tearing 
away  the  dull  grains  to  expose  fresh  cutting  edges  or  clear  away  the 
imbedded  material,  is  shown  in  Figure  2-5.  Before  the  teeth  are  cut, 
the  wheel  cutter  is  predressed  to  shape  the  cutter.  The  relation 
between  rotation  of  the  cradle  and  rotation  of  the  work  is  controlled 
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by  the  rotation  of  a  cam  which  rotates  about  a  fixed  center. 

The  formulae  derived  here  are  based  on  the  kinematics  of  two 
contacting  bodies.  There  is  no  relative  velocity  in  the  normal  direction 
at  the  cutting  point  and  the  teeth  are  assumed  to  be  rigid.  The 
formulae  are  written  in  vector  notation.  Most  of  these  vectors  are 
self-explanatory  except  for  the  following  two  terms: 

—  —  R  _ 

A(g,0)  ,  a  rotation  of  vector  A  through  an  angle  0  about  the 
unit  vector  g. 

- T 

(-,t,n)  ,  a  3  x  3  matrix  to  transform  a  vector  into  a  new 

coordinate  system  having  the  base  vectors  indicated 
in  the  parentheses.  Since  the  right-  hand  cartesian 
coordinate  is  used,  it  is  only  necessary  to  indicate 
two  base  vectors. 


2.2.1  Basic  Formulae 

To  describe  the  manufacturing  process,  one  needs  to  know  the 
relationship  among  the  rotating  elements,  and  the  relative  movement 
between  the  work  and  these  elements. 


2.2.2  Data  Defining  the  Gear  Tooth  Surface 


The  gear  tooth  surfaces  are  controlled  by  the  machine  settings. 

The  symbols  of  data  required  are  listed  in  Table  I.  Sample  data  are 
given  in  Appendix  A.  These  values  are  usually  tabulated  on  a  Gleason 
Gear  Engineering  Standard  form. 

A  spiral  bevel  gear  and  pinion  are  always  of  opposite  hand.  The 
formulae  considered  here  apply  to  a  set  consisting  of  a  left-hand  pinion  shown 
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TABLE  I 

SYMBOLS  OF  DATA  FOR  MACHINE  SETTING 


Number  of  teeth 

Face  width 

Shaft  angle 

Outer  cone  distance 

Addendum 

Dedendum 

Pitch  angle 

Face  angle 

Root  angle 

Spiral  angle 

Nominal  wheel  diameter 

Machine  center  to  back 

Sliding  base 

Blank  offset 

Cam  setting 

Eccentric  angle 

Cradle  angle 

Cam  guide  angle 

Standard  cam  setting 

Cam  pitch  radius 


Gear 


N 


g 

*». 

E 

mg 

T 

g 


g 

T 

go 

nu 

g 


n 


ig 


Pinion 

Concave  Convex 


F 

E 


M 


P 

xbP 

E 

mp 


bp 

a 

P 

T  „ 
pO 

nu 


n . 


ip 


X 

P 

E 

mp 

T 

P 


nu 

P 


n 

IP 


Index  interval 
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TABLE  I  (continued) 


Gear 

Pinion 

Concave 

Convex 

Dresser  block  angle 

E*d 

Outside  pressure  angle 

^dg 

*dp 

^dp 

Side  dresser  radial 

rdg 

rdp 

rdp 

Outside  dresser  arm  length 

Ldog 

^dop 

Ldop 

Outside  diamond  setting 

X, 

dog 

x, 

dop 

XJ 

dop 

Inside  dresser  arm  length 

T 

""dig 

L,. 

dip 

Ldip 

Inside  diamond  setting 

Xdig 

XdiP 

XdiP 

Side  dresser  axial 

Zdg 

Zdp 

Zdp 

Side  dresser  offset 

Edg 

Edp 

Edp 

in  Figure  2-6  and  right-hand  gear  shown  in  Figure  2-7.  The  drive  side 
is  assumed  to  be  the  concave  side  of  the  pinion  running  with  the  convex 
side  of  the  gear.  The  opposite  combination  is  the  coast  side.  Since 
the  gear  tooth  data  define  the  side  under  consideration,  these  general 
formulae  apply  to  either  the  drive  side  or  the  coast  side. 

2.2.3  Mechanism  of  Cam 

The  relative  rotation  of  the  work  and  cradle  is  controlled  by 
the  index  head  and  cam.  The  index  head  is  used  to  transfer  accurately 
the  rotation  ratio  between  the  cam  and  workpiece.  By  describing  the 
motion  of  the  cam  with  a  specific  index  head  setting,  one  can  obtain 
the  rotations  of  the  work  and  cradle,  and  their  angular  speed  ratio. 

The  follower  of  the  heart-shaped  cam  is  fixed  in  the  cradle,  and 
the  cam  rotation  center  is  fixed.  Figure  2-8  shows  the  position  for 
a  standard  cam  setting.  As  the  cam  rotates  about  its  own  rotation 
center,  the  fixed  follower  forces  the  cradle  to  rotate  about  its  center. 
The  required  rotation  of  the  work  is  obtained  by  using  the  proper 
indexing.  The  relative  motion  can  be  readily  changed  by  using  a 
different  cam  geometry,  cam  setting,  cam  guide  angle  and  index  interval. 
The  symbols  used  in  the  motion  formulae  are: 

£  cam  rotation 

A  work  rotation 

cradle  rotation 
a  cam  guide  angle 

A nu  change  in  cam  setting 

n  cam  pitch  radius 
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T  cam  setting 

T  standard  cam  setting 

o 

For  a  standard  heart-shaped  cam,  the  mechanism  of  motion  of  the 

cam  and  cradle  at  center  of  roll  can  be  thought  of  as  a  pair  of  internal 

gears  which  has  a  pitch  radius  for  the  pinion  and  DC  +  n^  for  the 

gear  (see  Figure  2-9).  If  the  cam  rotates  an  angle  eq,  the  cradle  will 

rotate  an  angle  n  /(n  +  DC)  times  e  . 

u  u  o 

The  change  in  cam  setting  Any  is  defined  as  the  actual  cam 

setting  T  minus  the  standard  cam  setting  Tq  (see  Figure  2-10)  .  The 

dashed  lines  are  for  the  standard  cam  setting  for  a  rotation  e  .  For  a 

o 

change  in  Anu»  the  following  relations  are  derived. 


ro  DC  +  n  ,  e  ° 
u  o 


An 


sin  *  sin  il<  +  -rrr  sin(e  -  ij>  ) 
o  DC  o  o 


e  »  e  +  ip  — 

o  o 


(2.1) 


(2.2) 


(2.3) 


where  ijr  and  e  are  the  actual  rotations  for  the  current  setting. 

A  center  line  of  the  follower  is  parallel  to  and  is  in  the  middle 
of  the  follower  surfaces.  The  cam  guide  angle  is  defined  as  the  angle 
between  this  center  line  and  the  line  from  the  cradle  center  to  the 
cam  rotation  center  before  the  cam  rotates  (see  Figure  2-lla) .  Because 
the  follower  is  fixed  in  the  cradle,  the  distance  from  the  cradle  center 
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to  the  center  line  of  the  follower  is  constant  and  is  always  equal  to  DC 


sin  a. 


Figure  2-llb  shows  the  relative  positions  with  a  cam  guide  angle. 


Using  eq  as  the  standard  cam  rotation  angle,  one  obtains 


*0  “  DC  +n  e0 
u 


(2.4) 


sin  1^2  =  sin  +  sin  a 


(2.5) 


♦  *  ll>2  -  a 


(2.6) 


e2  =  eQ  +  *2  ■  *0 


(2.7) 


e  »  e2  -  a 


(2.8) 


If  both  the  cam  setting  and  the  cam  guide  angle  are  changed,  one 


can  find  the  following  relations  based  on  the  geometry  shown  in  Fig.  2-12 


for  a  standard  cam  rotation  e  : 

o 


it  »  -  £ 

v0  DC  +  n  0 
u 


(2.9) 


sin  \\>2  =»  sin  +  sin  a  +  sin(eQ 


(2.10) 


4/  m  <l>2  ~  a 


(2.11) 


e2  "  eQ  +  *2  '  *2 


(2.12) 


e  -  e2  “  a 


(2.13) 


r.  r  v  V  v  v 
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A  =  —  e 
n 


The  ratio  between  the  cradle  and  the  work  rotation  becomes 


N. 


(2.14) 


d£  „  _n_  1M 
dX  "  nl  °M 


(2.15) 


where 


N. 


M 


n  cos  tj'n  +  An  cos(e.  - 
u  u  u  □ 


V 


DM  =  Nji  +  DC  cos  ij>2 


(2.16) 


(2.17) 


These  ratios  of  the  rotational  velocity  and  the  rotation  of  the  work  and  the 
cradle  are  necessary  for  computing  the  velocity  and  the  normal  direction 
at  the  contact  point.  Their  derivations  are  shown  in  Appendix  B. 

2.2.4  Basic  Machine  Set-up  and  Vector 
Coordinate  System 

The  cradle  settings  in  the  spiral  bevel  gear  cutting  machines 
of  the  Gleason  Works  are  shown  in  Fig.  2-13.  The  settings  are  fixed  during  the 
cutting  process.  From  this  figure,  the  following  relations  are  obtained. 


8. 


*P0 


16 

sin  -£■ 

(2.18) 

6 

_E 

2 

-  Q  +  360° 

P 

(2.19) 

6 

16 

sin  -£■ 

2 

(2.20) 

B 

qs 

_  _£ 

2 

(2.21) 
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Vectors  in  these  formulae  have  been  taken  with  respect  to  a 
set  of  right-hand  rectangular  cartesian  coordinates,  x,  y,  z  fixed  in 
the  machine  plane  viewed  from  the  front  of  the  cutter  toward  the  gear 
as  shown  in  Fig.  2-14. 

In  order  to  be  consistent  with  the  vector  formulae,  the  follow¬ 
ing  matrix  for  unit  vectors,  d,  e,  and  c  is  used. 

d  =  (1,  0,  0) 

7  =  (0,  1,  0)  (2.22) 

c  -  (0,  0,  1) 

2.2.5  Point  on  the  Gear 

In  cutting  the  gear,  the  moving  parts  include  the  rotating  work, 
cradle,  wheel  cutter,  cam,  and  dresser.  The  work  rotation,  cradle 
rotation  and  cam  rotation  are  related  to  each  other.  It  is  desirable 
to  choose  cam  rotation  as  the  independent  variable.  Thus,  there  are 
only  three  independent  rotations  which  define  a  specific  point  on  the 
gear  tooth  surface,  i.e.  the  cam,  dresser,  and  wheel  cutter. 

The  procedure  to  determine  a  point  on  the  gear  tooth  surface  is 
to  choose  arbitrary  values  for  the  cam  rotation  and  wheel  cutter 
rotation,  then  find  a  value  of  dresser  rotation  such  that  the  kinematic 
requirement  of  contact  is  satisfied,  i.e.,  the  relative  velocity  at  the 
cutting  point  is  perpendicular  to  the  normal  of  the  point  being  cut  in 
the  gear  tooth. 

Symbols  shown  in  Fig.  2-14'  for  all  position  vectors  are 
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defined  as: 

A  :  vector  from  the  machine  center  to  the  point  on  the  surface 

S  :  vector  from  the  machine  center  to  the  intersection  of  plane 

containing  the  cutter  tips  and  the  wheel  cutter  axis. 

R  :  vector  from  the  crossing  point  to  the  point  on  the  surface. 

B  :  vector  from  wheel  center  to  the  point  on  the  surface. 

r  :  vector  from  the  intersection  of  plane  containing  the  cutter  tips 

and  the  wheel  cutter  axis  to  cutter  tips. 

n  :  unit  normal  vector  to  the  cutter  blade. 

t  :  unit  vector  along  the  direction  of  dressing  diamond  movement, 

g  :  unit  vector  along  the  gear  axis, 

p  :  unit  vector  along  the  pinion  axis. 

The  initial  position  of  these  vectors  when  q  »  0,  0=0  are 

©  o 

g  =  (-cos  rr,  0,  -sin  r. ) 

Y  =  -y,  e 

'go  dg 

S  -  S  d 
go  g 

n  =  (0,  cos  $  ,  sin  $  ) 
go  g  g 

t  =  (0,  sin  4>  ,  -cos  $  ) 
go  V  g 

Because  the  convex  side  of  the  gear  tooth  is  formed  by  an  inside 
blade  of  the  cutter,  by  using  the  sign  convention  given  below,  the  same 
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formulae  (used  in  computer  programs)  can  be  used  for  both  sides: 


Drive  Side 


t>  ■  £$  -  b 

g  dg  dg 


X.  -  X,. 
dg  dig 


Ldg  Ldig 


Coast  Side 


♦»  '-♦dg 


X.  -  -  X, 
dg  dog 


Ldg  Ldog 


Reasonable  initial  values  of  cam  rotation  (standard  cam 

Og 

setting)  and  wheel  rotation  0^  are 


tn  -  0 

Og 


9g  qg0 


(2.23) 


(2.24) 


where  is  the  spiral  angle,  using  tn  from  cam  formulae  to  obtain 
M  ug 

and  ratio  of  cradle  rotation  to  work  rotation.  The  new  position 

g  8 

vectors  become 


q  *  q  +  ip 

4g  4g0  g 


S  -  S  _(c,  q  ) 
g  go"  4g 


r  ■  r  (c,  q  -9  ) 
g  go  4g  g' 


'i,  ■  Cg0^c*  VV 


;  ■  r  r  'tv.  r  i  ;vr 
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n,  *  n  (c,  q  -9  )R 
lg  g°  ’  4g  g 


R 


dlg  ‘  VC>  VV 


B  .  -  r  -  (0.2  -  a .  )  c  +  X.  n.  -  (2.0  -  E,  )  d. 
g0  g  dg  dg  lg  dg'  lg 


AgO  "  Sg  +  Bg0 

V  "  xg  1  ■ ¥mg  “  ‘  V  7  +  Ag0 


(2.25) 


where 


vector  from  wheel  center  to  dresser  center. 


vector  from  machine  center  to  dresser  center, 
vector  from  crossing  point  to  dresser  center. 


After  setting  certain  values  for  the  two  independent  variables 

and  9  ,  one  can  determine  the  dresser  angle  9  ,  that  satisfies  the 
0g  g  °  dg 

kinematic  requirement  of  contact  between  tooth  surface  and 

To  find  this  dresser  angle  9 ,  ,  an  initial  value  (9,  »  0)  is  fi- 

dg  dg 

assumed.  From  Fig.  2-14,  the  new  position  vectors  are 


d0  -  d,  (n  ,9.  ) 
2g  lg  lg  ag 


t2G  “  tlg{nlg’  9dg) 


A  *  A  n  +  L.  d, 
g  g0  dg  2g 


Rg  *  Rg0  +  Ldgd2g 


R 


(2.26) 
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B  =  B  n  +  L  d_ 
g  gO  dg  2g 


The  normal  of  the  cutter  blade  at  this  i  stant  (given  *  »  9  values)  is 

Og  eg 


n  =*{cxB)xt  } 
g  g  2g  unit 


(2.27) 


The  velocity  of  the  cutt  er  relative  to  the  work  at  the  point  being  cut 


-  m  —  —  _  — 

V  *  — —  cxA+gxR 

g  (cAj  g  6  g 


(2.28) 


The  product  n  •  V  is  checked  to  see  whether  or  not  it  is  zero;  if  not, 
g 

the  secant  method  is  used  to  iterate  0,  until  n  •  V  =  0  is  satisfied 

dg  g  g 

for  the  prescribed  values  of  e_  and  6  . 

Og  g 

The  tooth  surface  can  be  described  completely  with  the  formulae 
derived  so  far  by  varying  the  values  of  and  0^  independently. 

2.2.6  Point  on  the  Pinion 

The  procedure  to  find  the  pinion  tooth  surface  is  exactly  the  same  as 
that  for  the  gear  tooth  surface,  except  for  some  initial  vectors  which 
are  different.  From  Fig.  2-15,  the  mathematical  forms  are  described 
briefly  as  follows: 


Drive  Side 


r  r  , 

P  dp 


X,  * 
dp  dop 


L.  M  L. 
dp  dop 


Coast  Side 


♦p  '  -  r*dP 


X,  -  >  X.. 
dp  dip 


La  = 
dp  dip 


1  V' 
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The  initial  position  vectors  are: 

p  =  (-cos  Yr>  0,  -sin  y  ) 


S  -  S  d 
p0  p 


r  -  r,  e  (2.29) 

pO  dp 


n  *  (0,  -cos  <j>  ,  -sin  4  ) 
pO  P  p 


t  _  ■  (0,  -sin  $  ,  cos  ^  ) 
pO  P  P 


Assume  -  0,  8  =  q  _  -  one  can  calculate  X  ,4>  and  ratio  of 

Op  p  pO  M  P  P 

cradle  rotation  to  work  rotation  by  cam  formulae  described  earlier. 


qn  " 

p  pO  p 


s  =  s(c,  -q  ) 

p  pO  p 


r  *  r  (c,  8  -q  ) 
p  pO  P  P 


t.  =  t  (c,  8  -q  ) 
lp  pO  p  p 


n,  *  n  (c,  8  -q  ) 
lp  pO  P  P 


d  *  d(c,  8  -q  ) 
lp  P  P 


B  *  r  -  (0.2  -  ZJ  )  c  -  n  -  (2.0  -  E.  )  d, 
po  p  dp  dp  lp  dp  Ip 


>  ■  sp  +  Bp0 


R  x  p  +  E  e  -  X.  c  +  A 
pO  p  nip  bp 


(2.30) 
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Assume  a  value  of  dresser  angle  9 


dp 


d,  =  d.  (n  ,  e  y 

2p  lp  lp  dp 


C2p  =  ^p^lp’  9dpJ 


R 


A  =  A  +  L  & 
p  pO  dp  2p 


R  -R  +L.  d. 
p  pO  dp  2p 


(2.31) 


B  *  B  +  L .  d- 
p  pO  dp  2p 


n  *  {t  x  (c  x  B  )>  . 

p  2p  p  unit 


dt 

dA 


c  x  A  +  p  x  R 
P  P 


change  0^  until  V  •  n^  *  0. 

The  tooth  surface  of  pinion  is  thus  determined  by  changing  the 
values  of  and  9^  independently. 
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2.3  Tooth  Contact  Behavior  and  Motion  Graph 

The  only  tooth  curvature  of  involute  gears  with  an  infinite 
lengthwise  radius  of  curvature  is  the  tooth  profile  curvature,  so  the 
instantaneous  tooth  contact  pattern  will  be  a  straight  line  extending 
from  one  end  of  the  tooth  to  the  other.  If  elastic  deformation  is 
included,  the  contact  pattern  will  be  an  elongated  rectangle,  from  one 
end  of  the  tooth  to  the  other.  If  the  teeth  are  crowned,  lengthwise 
curvature  is  introduced.  This  results  in  a  point  contact  for  rigid 
teeth,  or  an  elliptical  contact  for  elastically  deformed  teeth. 

Theoretical  conjugate  bevel  gear  teeth  can  be  manufactured,  but 
they  are  sensitive  to  shaft  misalignment  and  housing  deformation.  Due  to 
this,  bevel  gears  are  produced  with  a  lengthwise  and  a  profilewise  mis¬ 
match.  Therefore,  the  instantaneous  tooth  contact  pattern  will  be  a 
theoretical  point  for  rigid  teeth  or  an  ellipse  when  elastic  deformation 
is  included. 

By  choosing  the  proper  amount  of  mismatch  in  both  directions  for 
a  particular  case,  edge  contact  and  excessive  contact  pressure  can  be 
avoided  in  the  presence  of  misalignment  due  to  shaft  and  housing 
deformations. 

The  purpose  of  the  motion  graph  is  to  describe  the  nature  of  the 
motion  transmitted  and  the  effects  of  contact  between  adjacent  pairs  of 
teeth  on  the  tooth  contact  pattern. 

The  gear  contact  problem  is  essentially  a  Hertzian  contact  prob¬ 
lem,  with  the  pressure  distribution  in  the  Hertzian  ellipse  a  main  concern 
in  EHD.  The  shape  of  the  Hertzian  ellipse  can  be  readily  determined 
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from  the  principal  radii  which  can  be  computed  analytically  from  the 
gear  teeth  geometry. 

2.3.1  Matching  a  Contact  Point 

A  method  to  find  the  contact  point  by  machine  setting  and  running 
position  is  derived  in  this  section. 

In  addition  to  the  rectangular  cartesian  coordinate  used  in 
Chapter  II,  a  point  on  the  tooth  surface  with  respect  to  the  tooth  axis 
and  the  crossing  point  can  be  specified  in  the  axial  plane  by: 

L  :  the  distance  along  the  tooth  axis  from  the  crossing  point. 

R  :  the  radius  from  the  tooth  axis. 

The  primary  effect  of  misalignment  is  to  shift  the  tooth  contact 
to  the  boundaries  of  the  tooth  surface,  resulting  in  edge  contact  and 
extreme  loading  conditions.  Figure  2-16  shows  a  pair  of 
spiral  bevel  gears;  the  following  errors  in  relative  running  position  of 
the  pinion  with  respect  to  gear  are  discussed: 

:  running  shaft  angle. 

X  :  pinion  apex  withdraw -which  is  the  linear  displacement  in 
pr 

the  H-direction  (along  the  pinion  axis). 

X  :  gear  apex  withdraw -which  is  the  linear  displacement  in  the 

gr 

Z-direction  (along  the  gear  axis) . 

Er  :  offset-which  is  the  linear  displacement  in  the  V-direction. 


Consider  a  point  on  the  gear  tooth  surface  obtained  by  assumed 

values  of  e  and  0  as  described  in  Section  2.2,  in  order  to  match 
Og  g 

the  gear  axis  from  its  manufacturing  position  to  a  proper  running 
position  relative  to  the  pinion  axis.  Suppose  one  lets  the  manufacturing 
axis  of  the  pinion  coincide  with  the  pinion  axis  at  running  position. 

Then  the  gear  axis  at  running  position  is  as  shown  in  Figure  2-17,  and 
the  new  position  vectors  relative  to  the  new  gear  axis  are: 


u  =  x  +  n  +  7T 
g  g 


az  -  z  +  r  +  t 

r  r  r 


•4> 

gl 


gl 


gl 


'gl 


g 


'2g 


.  R  —  .R 
(g,  y)  (e,  AZ ) 


(2.32) 


(2.33) 


where  y  is  the  desired  rotation  angle  of  the  gear  from  a  specified 
g 

reference.  In  the  computer  program,  n  is  set  at  zero  to  find  the  first 

g 

contact  point,  and  then  is  increased  to  obtain  the  successive  contact 
points . 


Assuming  that  a  point  P  on  the  gear  surface  obtained  by  a  pair 
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of  values  e  and  8  is  a  contact  point  under  running  conditions,  the 
Og  g 

position  vector  for  this  point  expressed  in  the  pinion  axis  can  be 
obtained  by  placing  the  gear  in  running  position  relative  to  the  pinion 
as  shown  in  Figure  2-18 


R  =  X  p  +  E  e  -  X  g  +R 
pi  pr  r  gr  1  gl 


(2.34) 


One  can  define  this  contact  point  P  in  the  pinion  without  con¬ 
sidering  its  angular  position  by 


L  =  -R  •  p 
pi  pi 


V  'V  x  p| 


(2.35) 


In  the  manufacturing  procedure,  there  is  a  set  of  particular 

values  of  and  6^  to  produce  a  point  in  the  pinion  with  the  same 

values  of  L  and  R  as  those  of  L  ,  and  R  , . 

P  P  pi  pi 

The  mathematical  way  to  find  these  particular  values  of  e  and 

9  is  described  as  follows. 

P 

Assuming  arbitrarily  and  0^  and  using  the  method  described 
in  section  2.2,  one  can  get  R^  which  is  the  vector  from  crossing  point 
to  the  contact  point  Q  between  cutter  blade  and  work.  The  vector  R 


can  be  specified  by 


R  =  |R  x  p 
P  P 


(2.36) 
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Iterate  e„  and  0  simultaneously  until  L  ■  L  ,  and  R  ■  R  ,  that  is. 
Op  P  p  pi  p  pi 


the  same  distance  around  the  pinion  axis  from  crossing  point  is  obtained 
for  point  Q  and  point  P. 

Then  rotate  the  pinion  axis  to  make  point  Q  coincide  with  the 


contact  point  P;  the  rotation  angle  required.  An  ,  is  shown  in  Figure 


2-19. 


An, 


sin 


(R  ,  x  p) 
_ 21 _ 


(2.37) 


Point  P  is  obtained  under  the  assumption  that  arbitrary  values 


of  en  and  0  were  made.  Under  a  given  tooth  profile  and  running  posi- 
8  8 


tion,  not  every  point  in  the  tooth  could  be  a  common  contact  point.  A 
common  contact  point  must  have  the  same  normal  direction  for  both  pinion 
and  gear. 

Assuming  that  a  point  P  is  the  common  contact  point,  the  normal 


direction  at  this  point  on  the  gear  for  and  0  is  n  , ,  on  the 

Og  g  gl’ 


(2.38) 


where  n  is  obtained  at  point  Q.  n  ,  and  n  ,  are  compared  to  see  if 
P  gl  Pi 


they  agree  to  within  a  specified  error.  If  not,  the  secant  method  is 


used  to  to  change  e  and  0  simultaneously  to  find  another  contact  point 

°8  8 


until  the  two  normal  directions  agree.  This  procedure  yields  a  point  of 
contact  between  gear  and  pinion. 
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In  fully  matched  gears,  the  constant  velocity  ratio  is  main¬ 
tained,  In  single  mismatch  case  (lengthwise  mismatch)  the  theoretical 
contact  point  is  a  point  of  the  original  surface,  so  the  uniform  velocity 
will  be  maintained.  There  is  only  one  theoretical  contact  point  in 
combined  mismatch,  the  point  of  the  original  surface;  therefore,  the 
velocity  ratio  will  not  remain  constant,  and  there  is  a  theoretical 
instantaneous  velocity  change  at  each  changeover. 

The  motion  graph  shows  the  error  in  angular  displacement  of  the  driven 
member  against  rotation  of  the  driver.  It  is  assumed  throughout  that 
the  pinion  is  the  driving  member  and  that  it  is  rotating  at  a  constant 
velocity. 

2.3.2  Condition  of  Contact 

The  surface  topography  of  spiral  bevel  gears  cannot  be  readily 
expressed  mathematically  by  an  analytical  function.  However,  one  can 
use  second-order  surface  approximations  to  describe  the  surface  around  a 
point  of  contact.  By  this  approximation  the  principal  radii  and  relative 
surface  velocity  at  the  contact  point  can  be  analytically  determined. 
Approximating  a  surface  to  the  second  order  gives 

z  =  0.5  a  x2  +  b  x  y  +  0.5  c  y2  (2.39) 

The  numeric  coefficients  were  chosen  so  that 


a 


curvature  in  x-direction 
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twisting 


(2.40) 


curvature  in  y-direction 


express  this  surface  equation  in  compact  form  (s)  = 


be 


Except  for  the  tooth  breakage,  the  tooth  failure  modes  are 


mainly  due  to  local  contact  pressure  which  causes  surface  pitting  or 


scoring  at  high  sliding  velocities.  Therefore,  to  develop  reliable 


formulae  for  gear  design  it  is  necessary  to  have  knowledge  of  contact 


pressure  and  surface  velocity  on  the  gear  tooth. 


2. 3. 2.1  Surface  Coefficients 


Figure  1-20  shows  the  surface  at  the  point  of  contact  for  the  gear 


The  tangent  plane  coordinate  system  is  chosen  in  such  a  way  that  the 


plane  containing  the  normal  direction  vector  of  contact  point  (n  )  and 


the  vector  c  of  wheel  axis  is  the  yz  plane  of  the  tangent  plane  coordinate 


system  as  shown.  Let  t  coincide  with  y  axis  and  n  with  z  axis,  so 
3  g  g 


the  vector  t  x  n  is  in  the  direction  of  the  x-axis. 
g  g 


In  this  chosen  tangent  plane  coordinate  system,  there  is  no 


twisting  for  the  surface  of  the  generator  which  is  the  wheel  surface 


around  the  point  of  contact  (b  =0)  and  the  curvature  in  the  x-axis 

swg 


direction  of  this  generator  surface  (a  )  can  be  derived  based  on  the 

SOOg 


geometry  shown  in  Figure  2-20. 


sin  \p  =  n  •  c 
wg  g 


(2.41) 


V  i  l 
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but 


hence 


r  =  B  x  c  sec  ip 
ng  g  1  <og 


stog  r 


(2.42) 


(2.43) 


ng 


2-21. 


The  curvature  in  the  y-axis  direction  (c  )  is  shown  in  Fig. 

siog  6 


t  =  n  tan  -  c  sec  <p 

g  g  log  (Og 


sin  i  -  n  •  (t„  x  n  ) 
g  lg  2g  g 


sin  a  *  t_  •  (t  x  n  ) 
wg  2g  g  g 


x  =  L.  sin  a  9  +  - 

dg  log 


y  =  L ,  cos  a  9  + - 

dg  tog 


A  =  0.5  L ,  (sin  i  )02  -  0.5x2  — 
dg  g  rng 


=  0.5  92[-L,  sin  i - ^  sin2  a  ] 

dg  g  r  iog 

ng 


(2.44) 


(2.45) 


(2.46) 


(2.47) 


A  =  0.5  y 


2  1 


iog 


1  1  Sln  1e  1.2 

-  =  - - —  tan  a 

iog  dg  cos  a  ng 

U)g 


cog 


(2.48) 
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1 

c  =  - 

S(Og  P 

OJg 


(2.49) 


where  p  is  the  radius  of  curvature  in  y-axis  direction. 

cog 

In  the  actual  running  position,  the  z-axis  of  the  tangent  plane 
is  n  ,  and  the  y-axis  is 

gl 


tgt  =  tg(g.V»)R  (e,  AZ)R 


(2.50) 


For  the  pinion,  the  surface  of  the  generator  around  the  point 

of  contact  is  shown  in  Fig.  2-22,  and  surface  coefficients  a  b 

scop  »  scop 

and  c  can  then  be  obtained: 
scop 
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a 

sup 


1 


r 

np 


b  -  0  (2.57) 

sup 


c 

sup 


1 


p 


up 


In  the  real  running  position,  the  tangent  plane  coordinate 
in  the  z-axis  direction  is 


and  in  the  y-axis  direction  is 


(2.58) 


-■  t  (p,  An  )R 
P  P 


(2.59) 


One  uses  the  vectorial  velocity  and  acceleration  of  a  known 
point  on  the  wheel  and  angular  velocity  vector  of  the  wheel  about  that 
point  to  describe  the  movement  of  the  generator.  This  point  is  chosen  as 
a  point  of  contact  at  a  given  instant.  The  motion  parameters  are: 


V  :  linear  velocity  vector. 

a  :  linear  acceleration  vector. 


u  :  angular  velocity  vector. 
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The  above  are  relative  motion  vectors  denoting  the  motion  of  a  moving 
generator  relative  to  a  generatee. 

The  values  of  these  motion  parameters  can  be  computed  for  the 
gear  as  follows  (see  Fig.  2-23  for  direction  of  vectors): 


.  Ot  = 

dV 

d\ 

O  c8 

o 

l«*J 

eg 


U)  =  c  +  g 
g  eg 


E  =  R  -  A 
g  g  g 


(2.60) 


(2.61) 


(2.62) 


V  =>  g  x  E  +o)  xA„  (2.63) 

g  g  2 

a  =■  g  x  (g  x  E  )  +  (w  I  x  c  +  a  c)  x  A  (2.64) 

o  ©  eg  Cg  o 

+  oj  x  (a:  x  A  )  (2.64) 

g  g  g 


Transforming  these  vectors  to  tangent  plane  coordinates: 


u)  *w(-,  t,n)^ 
tg  g  g  g 


vt8  ■  v*-  v  v1 


“«;  +  V'-  's'  V 


(2.65) 
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Similarly,  for  the  pinion,  one  obtains 


W  =*  u>  c  +  p 
P  cp 


E  -  R  -  A 
P  P  P 


V  =pxE  +o)  xA 
P  P  P  P 


a  *  p  x  (p  x  E  )  +  (u  pxc  +  a  c)xA 
P  P  cp  cp  p 


+  a)  x  (ui  x  A  ) 
P  P  P 


Transforming  to  the  tangent  plane  coordinate  system. 


a  -  a  (-,  t  ,  n  )T 
cp  P  P  P 


(2.66) 

(2.67) 

(2.68) 

(2.69) 

(2.70) 

(2.71) 


The  problem  is,  then,  to  determine  the  three  components  of  the 
difference  surface  matrix  (Aa,  Ab,  Ac)  in  terms  of  the  given  surface 
components  and  the  motion  parameters.  From  2.12,  the  following  relations  are 
obtained  for  the  gear  (subscripts  x,y,z  denote  the  component  of  a  vector): 
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a  +  V  a)  -  V  ui 
v  _ _ ££z _ tpx  tpy  tpy  tpx 

p  w  -  c  V 

tpx  sup  tpy 


u  -  c  V 

Ac _ £e* — gae. tpy _ 

p  v  +  V  -  V  tan  a 
p  tpy  tpx  sp 


Ab  =*  -Ac  tan  a 
p  p  sp 


Aa  =  -Ab  tan  a 
P  P  sp 


(2.79) 


(2.80) 


The  coefficients  of  the  surface  of  the  pinion  in  the  pinion  tangent  plane 
coordinate  system  are 


a 


sp 


a 

stop 


+  Aap 


b 

sp 


*  Ab 

P 


(2.81) 


c  *  c  +  Ac 
sp  sup  p 


It  is  necessary  to  rotate  the  pinion  tangent  plane  coordinate  system  to 
coincide  with  the  gear  tangent  plane  coordinate  system  in  order  to  express 
the  coefficients  of  both  surfaces  in  the  same  coordinate  system.  The  angle 
required  for  rotation  is  A  which  can  be  found  from 

sin  A  =  tpt  •  (tgt  x  ngl)  (2.82) 


and  then  to  transform  the  coefficients  of  the  pinion  surface  into  this 
coordinate  system, 
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2  2 
a  *  a  cos  A  +  2  b  cos  A  sin  A  +  c  sin  A 
spg  sp  sp  sp 


bspg  =  (csp  ~  asp)  cos  A  sin  A  +  bsp(cos2  A  "  sin2  A>  (2.83) 


2  2 
c  _  -  c  cos  A  -  2  b„„  cos  A  sin  A  +  a  sin  A 
spg  sp  sp  sp 


2. 3. 2. 2  Geometric  Contact  Ellipse 

Using  a  value  of  0.00025"  for  the  separation  at  the  extremities 
of  the  contact  and  the  surface  curvatures  of  each  member  at  the  point  of 
contact,  the  size  and  direction  of  the  contact  ellipse  are  determined 
as  (Fig.  2-24). 


Aa  *  a  -  a 

spg  sg 


Ab  =  b  -  b 

spg  sg 


AC  =  c _ -  c 

spg  sg 


(2.84) 


tan  2g 


2Ab 


Aa  -  Ac 


The  two  values  of  g  are  the  principal  directions.  Assuming  that  %2  *-s 
first  or  fourth  quadrant,  the  curvatures  in  the  g^  and  g^  directions  are: 
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1  £ 

—  =  Aa  cos  g1  +  2  Ab  cos  g  sin  g.  + 

p  ^  L  X  X 


A  '  2 

Ac  sin  g 


(2.85) 


1  2  2 
—  =  Aa  cos  g^  +  2  Ab  cos  g^  sin  §2  +  Ac  sin  §2 


and  p^  are  the  principal  radii  of  the  difference  surface  of  gear  and 
pinion  surface,  the  lengths  of  the  semi-axes  of  the  contact  ellipse  in  the 
tangent  plane  are: 


V  0.0005  -py 


pl 

VO.  0005  — 

_ 

2  _L 

p2 


(2.86) 


2. 3. 2. 3  Velocity  of  Contact  Point 

Consider  that  the  gear  surface  is  stationary.  During  a  small  unit 
time,  the  contact  point  on  pinion  moves  from  P  to  P'  with  velocity  Vg 
due  to  rotation  of  the  pinion.  During  this  same  time  interval,  the 
point  of  contact  moves  on  the  pinion  surface  to  P"  with  a  velocity  v 
(Fig.  2-25),  the  velocity  of  contact  point  is  the  sum  of  these  two 
parts . 


u  =  V  +  v 
s 


(2.87) 
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The  pinion  point  velocity  relative  to  the  gear  for  unit  pinion 


rotational  velocity  is 


Vsl  =  mi(Rgl  X  g)  +  (Rpl  x  p) 


Transforming  this  vector  into  the  common  tangent  plane 


V  =  V  (-,  t  .  n  .)T 
s  si  gt’  gl' 


The  angular  velocity  of  the  pinion  body  relative  to  the  gear  is 


=  -p  -  m±  g 


Transforming  this  vector  into  the  tangent  plane, 


“  ‘  V*-  V  V1 


(2.88) 


(2.89) 


(2.90) 


(2.91) 


The  surface  coefficients  of  gear  and  pinion  in  this  tangent  plane  are 


(s  )  = 

g 


a  b 
sg  sg 


b  c 
sg  sg 


a  b 
spg  spg 


b  c 
spg  spg 


(2.92) 


Since  p"  is  to  be  a  point  of  contact,  the  rate  of  change  of 
the  normals  of  the  gear  surface  and  the  pinion  must  be  the  same: 


:  *  •  |->A  •  «  v  V  V  V  v  v- 
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dn  — .  \ 

—  =  “  u(s) 
dt  g 


=  -  v(s  )  +  03  x  n 
P 


for  the  gear 


for  the  pinion 


where  03  x  n  is  due  to  the  rotation  of  the  pinion  surface,  and  -v(s  )  is  due 

P 

to  curvature  of  pinion  surface.  Hence, 


-u(s  )  =-(u  -  V  )  (s  )  +  03  x  n 
g  b  P 


u[(s  )  -  (s  )  ]  =  oj  x  n  +  V  (s  ) 
g  P  s  p 


u(As)  =  oi  x  n  +  V  (s  ) 
s  p 


u  =  [oixn  +  V  (s  )]  (As) 
s  p 


(2.93) 


V  T  t' 


s’* 
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2.4 _ Results  and  Discussion 

One  typical  set  of  spiral  gears  currently  being  tested  at  NASA 
Lewis  Research  Center  was  used  to  illustrate  the  numerical  calculations. 

The  machine  settings,  which  are  provided  by  Gleason  Works,  are  shown  in 
Appendix  A.  The  points  on  the  tooth  surface  were  obtained  by  procedures 
described  in  Section  2.2.  The  contact  condition  was  determined  by 
the  method  shown  in  Section  2.3.  A  series  of  solutions  were  obtained 
to  simulate  the  spiral  bevel  gears  subjected  to  a  range  of  operating 
conditions,  and  results  are  presented  here.  In  order  to  obtain  a 
unified  view,  all  the  curves  are  plotted  in  the  axial  plane. 

2.4.1  Tooth  Profile 

The  methods  in  Section  2.2,  which  use  actual  cutting  motions 
of  the  gear  cutting  tools,  result  in  an  exact  description  of  the  tooth- 
forms  as  functions  of  spiral  angle,  blade-edge  radiu.  ,  etc. 

A  typical  tooth  profile  for  the  gear  shown  in  Fig.  2-26  was 
plotted  in  Fig.  2-27.  For  the  purpose  of  applying  FEM  in  gear  design, 
a  typical  finite  element  model  for  a  single  gear  tooth  was  shown  in 
Fig.  2-28. 

2.4.2  Contact  Pattern 

The  spiral  bevel  gears  are  relieved  both  profilewise  and  length¬ 
wise,  so  that  there  is  only  a  single  contact  point  in  the  center  of  the  tooth. 
This  combined  mismatch  permits  the  gear  pair  to  tolerate  displacements 
under  load  and  assembly  errors.  Theoretical  point  contact  is,  of  course, 
only  an  imaginary  concept.  Under  even  light  loads,  tooth  surfaces 
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compress  and  deflect  to  cause  surface  contact.  This  surface  contact 
area  is  roughly  elliptical.  The  sum  of  those  elliptical  contact  areas 
along  the  path  of  contact  is  called  "Tooth  Bearing."  Correctly  made 
teeth  will  show  a  tooth  bearing  centrally  located  on  the  tooth  surface, 
when  the  gears  are  in  their  standard  relative  running  positions.  The 
amount  of  lengthwise  mismatch  is  readily  observed  by  the  length  of  the 
tooth  bearing.  Profile  mismatch  is  not  observed  as  readily  owing  to 
the  fact  that  the  tooth  height  is  relatively  small;  nevertheless,  the 
existence  of  profile  mismatch  is  evidenced  by  the  lack  of  any  heavy 
concentration  of  bearing  at  the  tips  of  the  teeth. 

The  results  illustrated  here  were  obtained  by  analyzing  the 

effects  of  running  position  errors  on  the  path  of  contact.  The  errors 

discussed  here  include  pinion  apex  withdraw  (X^) >  gear  apex  withdraw 

(X^) ,  running  shaft  angle  (Ef)  and  running  offset  (Er) .  Their  effects 

are  calculated  by  changing  the  value  of  9  single  term  while  keeping  the 

others  unchanged,  and  these  results  are  shown  from  Fig.  2-29  to  Fig. 

2-33.  The  path  of  contact  moves  to  the  toe  as  the  pinion  apex  withdraw 

increases  and  moves  to  the  heel  as  it  decreases.  The  shape  of  the  path  of 

contact  remains  the  same,  but  is  inclined  more  as  it  approaches  both 

sides.  The  gear  apex  withdraw  will  produce  the  same  effect,  but  in  the 

opposite  direction.  As  X^r  increases,  the  path  of  contact  moves  to  the 

heel,  and  the  opposite  occurs  as  X  decreases.  Both  of  these 

gr 

have  as  their  chief  effect  a  pressure  angle  error;  the  tooth  contact 
will  move  toward  the  tooth  tips  on  the  gear  tooth  and  toward  the  tooth 
flanks  on  the  pinion.  When  the  shaft  angle  becomes  larger,  the  contact  on 
both  the  gear  and  the  pinion  shifts  to  the  toe.  At  the  same  time,  the  contact 
moves  toward  the  tooth  tips.  The  contact  goes  the  opposite  way  when  the 


shaft  angle  decreases.  This  effect  is  expected.  The  principal  effect 
of  the  offset  is  to  tilt  the  tooth  in  its  mating  slot,  so  that  the  con¬ 
tact  shifts  to  the  toe  on  one  side  and  to  the  heel  on  the  other  side. 

Usually,  the  running  position  errors  are  introduced  simul¬ 
taneously;  the  results  on  the  drive  side  are  investigated  for  the  sliding 
velocity  and  radii  of  the  difference  surface  at  the  contact  point  with 
three  separate  sets  of  errors  which  result  in  the  center  contact,  heel 
contact,  and  toe  contact.  The  tooth  bearing  for  these  three  contact 
cases  is  plotted  respectively  on  the  same  graph.  Assuming  the  combined 
Hertzian  deformation  to  be  0.00025", which  corresponds  approximately  to 
to  the  case  of  a  pair  of  lightly  loaded  gears,  and^using  the  curvatures 
of  the  contacting  surfaces,  the  size  and  direction  of  the  contact 
ellipse  are  determined.  The  larger  the  radius,  the  bigger  the  ellipse. 
The  sliding  velocity  has  also  been  plotted  at  selected  points.  See 
Fig.  2-34  and  Fig.  2-35. 

2.4.3  Motion  Graphs 

The  rotation  errors  of  the  driven  member  against  the  driving 
member  for  the  entire  duration  of  contact  for  several  adjacent  meshing 
pairs  were  plotted.  The  corresponding  velocity  errors  were  shown  in 
the  same  graph.  Fig.  2-36  shows  the  curves  on  the  drive  side  with  a 
running  position  which  results  in  a  center  contact.  Fig.  2-37  shows 
the  results  on  the  coast  side  with  the  same  running  position.  Figs. 


2-38  and  2-39  illustrate  the  contact  patterns  for  the  motion  curves 
shown  in  Figs.  2-36  and  2-37. 
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The  velocity  curves  in  those  figures  do  not  intersect;  at  the 
transfer  points  there  is  a  sudden  jump  in  the  velocity  with  which  the 
driving  tooth  strikes  its  mate.  This  value  is,  therefore,  a  measure 
of  the  impact  occurring  at  each  tooth  mesh.  As  mismatch  is  increased, 
the  impact  increases  in  magnitude,  resulting  in  noise  and  stress.  It 
is  therefore  essential  that  mismatch  be  kept  to  the  smallest  value 
compatible  with  the  adjustability  to  the  running  position  errors.  This 
concept  of  motion  errors  is  only  applied  to  the  case  under  a  very  light 
load.  The  elastic  deformation  due  to  moderate  load  will  smooth  out  these 
errors. 
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A.XlAI.  SECTION  \l 


Fig.  2-2.  Sectional  View  in  The  Axial  Plane  For  A  Pair 
of  Spiral  Bevel  Gears 
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Fig.  2-3.  Sectional  View  in  The  Transverse  Plane  For  a  Pair 
of  Spiral  Bevel  Gears 
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Wheel 

Axis 


Fig.  2-5.  Movement  of  the  Dresser 


RIGHT  HAND 


LEFT  HAND 


Fig.  2- 


7.  Spiral  Bevel  Gears 


Cam  Rotation  Center 


DC 


Cradle  Center 


Fig.  2-9.  The  Movement  of  Standard  Setting.- 


Fig.  2-10.  The  Movement  of  Cam  With  a  Change  in  Cam  Setting 
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Fig.  2-12.  The  Movement  of  Cara  With  Change  in  Setting  ,aJ?a 
and  Cam  Guide  Angle 
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Pig.  2-17.  Relation  Vectors  Between  Running  Position  And 
Manufacturing  Position 


Fig.  2-18.  Position  Vectors  in  Running  Position  Relative  to 
Pinion  Axis 


Fig.  2-19,  The  Angle  Between  R 
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Dresser  \ 

Axis  ^ 

Fig.  2-20. 

Wheel  Surface  Curvature  (  gear  ) 

in  X  Direction. 

Wheel  Axis 


wheel  Center 
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on  Vectors (Considering  Gear  Stationary 


Fig.  2-25.  Surface  Velocity 
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Fig.  2-27.  Tooth  Profile  Plotted  by  Calculated  Points 
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CHAPTER  III 

DYNAMICS  OF  SPIRAL  BEVEL  GEARS 

3. 1 _ Introduction 

One  of  the  main  uncertainties  in  gear  failure  analysis  is  the 
load  imposed  on  the  gear  teeth  at  high  speeds  where  the  inertia  forces 
of  the  gear  wheel  and  shaft  become  significant.  The  dynamic  response  of 
the  operating  system  plays  an  important  role  in  gear  design  because  the 
movement  of  the  gear  body  will  shift  the  contact  bearing  to  an  undesirable 
position  as  a  result  of  the  misalignment  of  shaft.  The  dynamic  load  in 
gearing  has  received  continuous  attention  in  the  past,  and  many  studies 
have  been  reported  (1-9).  As  far  as  the  dynamic  response  of  bevel  gears 
is  concerned,  most  of  the  studies  have  been  conducted  experimentally  in  the 
laboratory.  One  of  the  recent  measurements  is  by  Texauchi  and  Fuji!  (10) 
in  which  the  dynamic  load  of  gear  teeth  of  straight  bevel  gears  is 
obtained  bv  using  a  power  circulating  gear  testing  machine.  There  is 
a  lack  of  analytical  study  of  the  vibration  of  spiral  bevel  gears.  For 
instance,  some  industrial  practice  still  depends  upon  Buckingham's  formulae 
or  upon  empirical  factors  for  designing  spiral  bevel  gears. 

Analytical  studies  of  spiral  bevel  gears  are  extreme Iv  complex 
because:  1.  trie  tooth  profile  of  spiral  bevel  gears  is  generated  bv  -tie  cutting 

machine ,  and  it  is  difficult  to  describe  the  tooth  sur: ace, which  is  me. ■ 
more  complicated  than  that  of  spur  and  he  1  i  ai  gears.  the  toot:-,  -nt  .  * 
patterns  -f  spiral  bevel  gears  van  with  the  assemble  errors  a  .sec.  bv 
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the  apex  error,  shaft  alignment,  etc.;  3.  being  different  from  spur 
gears,  the  spiral  bevel  gears  have  also  lateral  vibration  in  addition 
to  rotational  vibration.  In  this  report,  efforts  are  made  to  explore 
the  effects  of  all  these  factors  on  the  dynamic  behavior  of  the  spiral 
bevel  gears. 
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3. 2 _ Formulation 

The  main  objective  in  this  section  is  to  formulate  a  dynamic 
model  and  equations  of  motion  to  simulate,  in  the  steadv  state,  the  periodic 
motion  of  both  the  pinion  and  the  gear  as  well  as  the  tooth  loan  during 
a  tvpical  cvcle,  in  which  a  pair  of  teeth  traverses  through  the  zone  or 
action  from  point  A  to  point  C  as  shown  in  Fig.  3-1. 

In  developing  the  equations  of  motion,  the  pinion  anc  the  gear  were 
assumed  to  be  rigid  bodies  each  having  six  degrees  of  :reedom.  The 
supporting  bearing  was  assumed  flexible  with  gnown  spring  st  itfness  At 
the  contact  or  each  pair  c  f  teeth  m  the  zone  •  t  a.  t  ion .  r  he  teeth  were 
assumed  oe  onneoted  bv  a  linear  spring  wit:  -  -r lent  «»a  no'tti .  : 

the  ont  a  t  p  • ;  n  t  and  nas  a  st  t  f  t  ness  t  •  t*e  iet  ermine  :  «»tar  ate  .  i 

"  e  .  erne  n  :  Th  -  Je  1 

The  eqr.at  lorss  ’  »”  ;  -r'  ■*  ,  •  :  g  . -  •  ••  .  •  :v-«-  :  '.srr-  l’- 

see  rig 

w  - 

W*>  ,  •  •  1  '  .•»*•*'-(  ■  •  „  M 

•»  - 
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The  relieving  derivation  is  described  separately  for  pinion  and  gear. 


i . 2 . 1  Pinion 

Let  the  moving  coordinate  axes  xf,  y' ,  z'  be  the  principal  axes 

P  P  P 

the  inertia  of  the  body  so  that  the  axis  of  the  rotation  of  the  body 

will  airert  alone  one  of  these  axes,  which  is  named  the  dominant  axis. 

Tne  pinion  shaft  is  not  fixed  in  the  moving  coordinate,  but  rotates 

it>.  at  tne  dominant  axis  y'  at  an  average  velocity  io  which  is  the 

p  P 


perationai  rotation  velocity  of  the  shaft.  XYZ  axes  are  fixed  co- 
rzmates  in  space  (see  Fig.  3-1).  In  the  case  of  the  pinion,  the 
ruminant  axis  is  along  the  Y  axis,  hence  the  angular  momentum  of  the 
n  about  the  mass  center  is 


-  ,  i'  +  I  ,  (6  , 

x  p  p  y  p  y  p 


+  u)  )  j '  +  I  ,  9  .  k' 
P  p  z  p  z  p  p 


(3.2) 


rate  o: 


range  of  in  the  moving  coordinate 


•  ,  i’  +  I,  9  ,  j'  +  I  ,  9  ,  k*  (3.3) 

xpp  YPypp  z  p  z  p  p 


:  x  (I,  9  ,  i’  +  I  ,  (9  ,  +«)j'  +1  ,  9  ,  k') 

p  xpxpp  yp  yp  p  p  zpzpp 


•  .  i'  *  "  ,  j'  +  e  ,  k’)  (3.4) 

xpp  y  p  p  z  p  p 


-  I  ,  a  9  ,  +  ( I  ,  -  I  ,  )  9  ,  9  ,  ]  .  , 

x  p  x  p  vppzp  z  p  yp  ypzp  l 
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+[1,0,1  j' 

y  p  y  p  J  p 

(3.8) 

+  1,9,  +  I  ,  a)  9  ,  1  k' 

L 2  p  z  p  y  p  p  x  pj  p 


where  M  is  number  of  bearing  reaction  force,  N  is  number  of  contact 

force,  and  x'  ,  v'  ,  and  z'  are  the  components  of  the  absolute 
cp  ' cp  cp 

acceleration  of  the  pinion  mass  center  in  the  x’,  y’  and  z'  direction. 

} .  »  .  *  >e  a  r 

These  results  are  similar  to  the  equations  for  the  pinion,  except 
•■lit  trie  tonil'  ant  axis  in  the  case  of  the  gear  is  X  axis  and  the  operational 
•it:  •  .  i  tv  :s  ,  the  equations  of  motion  for  the  gear  are: 

-  S 

»  <*’  i*  ♦  v'  j*  +  z’  k* )  (3.9) 

*  >'*  *  eg  g  eg  g 

if  *  Output  Torque 

xi  *  l 

*  K 

1  J . 10; 
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Some  preliminary  work  must  be  completed  before  the  twelve 
equations  in  the  system  can  be  solved.  These  are:  1)  expressing  the 
force  terms  and  Fci  in  terms  of  the  displacement  unknowns  in  the 
fixed  coordinate;  2)  expressing  the  displacement  unknowns  in  the 
fixed  coordinate  by  the  displacement  unknowns  in  x’y'z'  system  in  which 
the  equations  are  solved. 

Because  the  contact  force  is  assumed  to  always  be  in  the  normal 
direction  of  the  contact  point  on  the  tooth  surface,  the  scalar  quantity  of  the 
contact  force  can  be  used  effectively  in  the  derivation  if  the  direction 
vector  of  the  contact  point  is  known.  For  simplification,  the  matrix 
notation  is  used,  and  some  basic  column  matrices  are  introduced: 


{D' }  =  (x' ,  y',  z't  0  ,  ,  8  ,  0  ,  ) 

g  g  g  g  x'g’  y  g  z'g 


{D'}  =  (x' ,  y',  z',  9  ,  ,  6  ,  ,  9  ,  ) 

J p  p *  ^ p *  P*  X'p>  y'p  2 f  p 


(3.11) 

(3.12) 


which  are  the  displacement  unknowns  in  the  moving  coordinate  system, 
while 


{  D  } 

g 


(x  ,y  ,z  ,0  ,8  ,0  ) 

g  V  g  Kg*  ys  z8 


(3.13) 


(D)  *  (x  ,y  ,z  ,0  9,„>9  )' 

p  p  'p  p  xp,  yp  zp 


(3.14) 


are  the  displacement  components  in  the  fixed  coordinate. 

From  the  configuration  of  the  gear  system, the  displacement 
components  at  the  contact  point  i  in  the  normal  direction  (G^)  can  be 
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expressed  in  terms  of  the  rigid  body  displacements  as: 


(G  ,  > 

ni  g 


[DG 


iJg 


{D} 


(3.15) 


{G  }  =  [DG. ]  {D}  (3.16) 

ni  p  i  p  P 

and  the  bearing  reaction  force  at  bearing  j  can  also  be  expressed 
in  terms  of  rigid  body  displacements  by: 

(F  }  =-[DK  ]  {D}  (3.17) 

rJ  g  i  g  g 


(F  . }  -  -  [DK. ]  (D>  (3.18) 

rJ  P  J  P  P 

where  i  and  j  are  from  one  to  the  number  of  contact  forces  and  number  of 
bearings, respectively .  The  displacements  at  contact  point  i  due  to 
elastic  deformation  along  the  direction  of  contact  normal  are: 


(6  .) 

nci  g 
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where  i,  j  are  from  one  to  the  total  number  of  contact  points.  The 
unit  outward  normal  vector  at  contact  point  i,  is  denoted  nci.  The 
positive  value  of  Fci  is  taken  so  that  contact  force  is  acting  along 
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Che  sane  direction  as  the  outward  pointing  normal  vector.  Since  only  the 

scaLar  quantity  of  the  contact  force  is  considered  in  the  equation,  the 

relation  : F  . 1  »  \F  . }  is  alwavs  true  at  each  contact  point. 

ci  g  ci  p  -  K 

The  total  displacements  in  the  normal  direction,  whi  h  include 

both  the  rigid  body  and  the  elastic  deformations,  will  be  the  same  for 

the  two  contacting  bodies.  This  relation  gives: 


n  ,  •  ({<$  }  +  {G  }  )  -  -  n  *  ({6.}  +  {G  }  ) 

cig  i  g  i  g  cip  x  p  1  p 


(3.21) 


{5  .}  +  (G  }  -  -  {5  .}  -  (G  .} 

nci  g  ni  g  nci  p  ni  p 
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[DC.,]  ( F  . }  +  (G  }  =  -  [DC  ]  (F  }  -  (G  } 
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and 


-1. 


[F  ,}  =  -  ( [DC  ]  +  [DC..]  )  ({G  .}  +  (G  .)  ) 
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(3.23) 


(F  ,}  -  -  ([DC.  1  +  [DC,.]  )_1 ( [DG  ]  {D}  +  [DG  ]  {D}  )  (3.24) 

cig  ji  g  ji  P  i  g  g  ipp 


which  is  the  equation  relating  the  contact  force  and  the  rigid  body 
displacement  unknowns.  The  displacements  in  different  coordinate 
systems  can  be  transferred  from  one  to  another  by  means  of  the  following 
relations  (10), 


i  i*  v*  r  ir  r  v  *■» V*  • 

i  l  l  i  .*  *> 
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(3.25) 


1  <  y 


(3.26) 


,u,  r.k 


[9  1  [8  ]  {  y 


(3.27) 


u 

where  [9  ]  is  the  transformation  matrix  for  rotation  in  x'y'z’  coordinate. 


[9  ]  is  the  transformation  matrix  for  rotation  in  xyz  coordinate. 


Incorporating  all  the  above  force  and  displacement  relations 


into  the  equations  of  motion  for  the  pinion  and  gear  and  then  rearranging 


the  matrices,  one  can  obtain  the  final  form  of  the  twelve  equations. 


+  [c] 


+  [k] 


(3.28) 


The  details  of  matrix  [m],  [c]  and  [k]  are  derived  in  Appendix  C. 


After  examining  the  equations  of  motion,  it  was  found  that  the 


fourth  equation  governing  the  rotational  motion  along  the  dominant  axis 


't  the  gear,  ana  the  eleventh  equation  governing  the  ri’tdt  l^nal  motion 


aiiinK  the  dominant  axis  of  the  pinion  .  an  he  reJm  ed  to  a  single  equation  . 

The  reduced  equation  is  similar  to  that  describing  the  vibratorv  system  of 

two  masses  connected  bv  a  spring  In  this  equation,  it  is  convenient  to  convert 

the  rotational  motion  into  a  motion  along  an  artificial  line  of  action. 

Like  the  line  of  action  in  spur  gears,  the  artificial  line  of  action  In 
spiral  bevel  gears  is  defined  as  the  effective  component  of  y  x  n  to 
produce  the  torque  along  the  dominant  axis,  where  i  is  a  position  vector 
from  the  mass  center  to  the  contact  point,  and  n  is  a  unit  normal 
vector  at  the  contact  point. 

The  4th  and  11th  equations  of  motion  are  rewritten  here  as 


I, a,  ■  (I  y  .  x  F  ),  +  (Output  Torque)  ,  (3.29) 

X  g  x’g  cgi  cgi  x'g  X  g 


I  ,  a  ,  •  (l  Y  .  x  F  .),  +  (Input  Torque)  ,  (3.30) 

y  p  y  p  cpi  cpi  y  p  y  p 


The  following  new  variables  are  introduced: 


0  -  RBG  •  9  ,  (3.31) 

g  x'g 

9  *  RBP  •  9  ,  (3.32) 

P  y  P 

Letting  FMA  denote  the  magnitude  of  the  contact  force  and  Tc  be  the  static 
force,  the  equations  become 
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(3.36) 


RBC  and  RBP  are  treated  as  the  radii  of  the  effective  base  circle  of  spiral 
bevel  gears.  Eqs.  (3.29)  and  (3.30)  become 
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Defining 


J 


xg 


and 


yp 


RBP2 


one  obtains 


J  0  =  FMA  -  T 

xg  g  c 


J  0 
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(3.39) 


(3.40) 


(3.41) 

(3.42) 


'.itroducing  an  equivalent  mass  J  and  a  relative  displacement  *3 


Pg 

(3.43) 


4  a  4  -  4 

Pg  P  g 


the  reduced  equation,  including  damping  terms,  becomes 


J6  +  C  3  +  FMA  «  T 

Pg  pg  c 


(3.44) 


(3.45) 


After  the  relationship  between  the  11  displacement  unknowns  and  the 
magnitude  of  the  contact  force  is  rewritten,  the  system  can  be  reduced  to 
11  second  order  differential  equations  with  state  dependent  coefficients. 
These  equations  were  solved  numerically  by  the  Runge-Kutta  method.  The 
unknown  initial  conditions  (11  displacements  and  11  velocities)  can  be 
iterated  by  taking  the  previous  calculated  values  at  the  end  of  each 
iteration.  The  criterion  of  the  iteration  is  to  force  the  equation  to 
satisfy  the  periodic  contact  condition.  The  static  displacements  and 
zero  velocities  are  taken  as  first  guessed  values. 

To  implement  the  above  equations,  the  following  state  coeffi¬ 
cients  are  needed: 

(1)  The  tooth  contact  position  as  a  function  of  the  relative 
rigid  body  displacements  of  the  two  shafts, 

(2)  The  direction  of  the  normal  vector  at  the  contact  point, 

(3)  The  combined  stiffness  of  the  teeth  at  the  contact  point. 
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The  teeth  contact  position  and  the  direction  of  the  normal  vector  at  the 
contact  point  can  be  obtained  from  Chapter  II.  Because  cf  the  geometric 

complexity  of  spiral  bevel  gears,  there  is  no  satisfactory 
formula  available  to  calculate  the  deformation  of  the  tooth  surface. 

A  finite  element  program,  to  be  described  in  Section  3.3,  was  used  in 
this  study  to  calculate  the  deformation  due  to  a  unit  load  applied  at  a 
given  contact  point  of  spiral  bevel  gears. 
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3 .  3 _ Tooth  Deflection 

For  most  gear  applications,  the  contact  ratio  is  greater  than 
one,  that  is,  there  will  be  more  than  one  pair  of  teeth  in  contact 
during  some  portion  of  the  whole  engagement.  When  the  load  is  shared 
by  two  contacts,  it  cannot  be  assumed  that  the  load  is  distributed 
equally  among  the  pairs  of  teeth  in  contact  because  this  is  a  statically 
indeterminate  case.  Therefore,  one  must  consider  the  tooth  deflection 
under  the  load  for  each  pair  in  order  to  determine  the  load  sharing 
characteristics  among  the  pairs. 

Because  of  the  complexity  of  the  spiral  bevel  gear  geometry, 
there  are  no  existing  simplified  methods  for  calculating  the  suitable 
tooth  deflection.  In  order  to  investigate  the  system  response,  the  gear 
shaft  must  also  be  included  in  the  calculation  of  tooth  deflection. 

Some  of  the  recent  applications  of  finite  elements  in  determining 
tooth  deflection  can  be  found  in  Refs.  9  and  12,  where  it  was  shown  that 
more  accurate  results  can  be  obtained  by  using  the  finite  element  method. 
Most  of  this  work  dealt  with  two  dimensional  problems  and  did  not 
include  the  whole  gear  body.  Fig.  3-4  shows  a  typical  8-node  solid 
element  grid  pattern  for  a  gear  and  a  pinion  with  three  adjacent  teeth 
attached  to  the  gear  wheel  and  shaft.  Fig.  3-5  shows  the  central 
tooth  and  its  attached  ring  element  of  the  gear.  Fig.  3-6  shows  parts 
of  the  gear  shaft  and  the  gear  wheel .  Fig.  3-7  shows  the  entire  ring 
element  with  three  adjacent  teeth  of  the  pinion.  Fig.  3-8  shows  the 
elements  of  the  pinion  shaft.  The  central  tooth  was  subject  to  load  to 
calculate  the  deflection.  There  are  941  nodes,  562  elements  for  the 
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a  load  P  applied  at  anv  grid  t>oint  in  the  tooth  surtaoe.  Ft  this 
.analysis,  the  MAKC-CDC  program  was  used,  the  boundaries  were  considered 
to  be  fixed  for  ail  the  points  connected  to  thrust  bearing  to  eliminate 
rigid  body  displacement ,  and  the  boundaries  nodes  connected  to  the 
radial  bearing  was  allowed  to  move  in  the  direction  ot  dominant  axis. 

The  stiffness  at  grid  point  i  is  defined  as 


KS  .  *  ~  (3.46) 

gl  V 

KS  .  •  (3.47) 

pl  4pl 

The  stiffness  (KSg,  KSp)  of  the  point  other  than  the  grid  point  in  the 
tooth  surface  can  be  calculated  by  the  interpolation  method.  The 
details  are  shown  in  Appendix  D.  The  combined  stiffness  at  contact 
point  is  expressed  as 

KS  •  KS 

KS  =  KS§  +  KS?  (3.48) 
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ect  damping  rat ios  and  .. ant  act  ratios. 

.1  Dynamic  Lead  Variation 

For  a  constant  input  torque,  the  load  on  the  contact  point  of 
the  two  meshing  teeth  along  the  path  of  contact  is  not  constant.  This 
Load  variation  is  mainly  caused  by  the  following  factors: 

1.  The  variation  of  stiffness  along  the  contact  path. 

2.  The  transition  from  a  single  pair  of  contacts  to  a  double 
and  from  double  to  single. 

3.  The  effective  radius  is  not  constant  along  the  contact  path. 
Fig.  3-9  shows  the  variation  of  effective  radius  of  pinion.  Fig.  3-10 
shows  the  variation  and  jump  of  stiffness  for  the  transition  of 
contacts . 

The  main  excitation  to  the  gear  system  comes  from  the  periodical 
change  in  teeth  stiffness  due  to  the  alternating  engagement  of  single 
and  double  pairs  of  teeth.  The  frequency  of  this  excitation  force, 
expressed  as  a  meshing  f requency, depends  on  the  operating  speed. 
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Therefore,  it  dominates  the  resulting  mode  of  vibration.  Fig.  3-11  to 
Fig.  3-14  show  dynamic  load  variation  at  four  different  speeds  in  the 
case  of  central  contact,  that  is ,  when  the  contact  path  is  located 
centrally  between  the  toe  and  the  heel  of  the  tooth. 

Since  there  are  eleven  degrees  of  freedom  in  the  system,  eleven 
resonating  frequencies  of  the  system  should  exist.  In  the  low  speed 
region  where  the  excitation  frequency  from  the  change  of  stiffness  is 
much  lower  than  all  resonating  frequencies,  the  dynamic  load  response 
along  the  path  of  contact  is  somewhat  like  static  load  superimposed  fcv 
an  oscillatory  load  due  to  the  resonating  frequency  of  the  system. 

When  the  speed  is  near  the  resonance  region  (Fig.  3.11).  the 
dynamic  load  response  becomes  very  severe  (Fig.  3-12  and  Fig.  3-13  . 

The  maximum  dynamic  load  is  much  higher  than  the  static  load,  which  . 
the  case  when  overloading  occurs.  Sometimes  the  oscillation  o:  dvnarr. 
load  will  make  meshing  teeth  separate  when  the  load  becomes  necat  :  e 
and  thus  generates  noise  and  surface  fatigue. 

As  the  speed  increases  beyond  the  zone  of  resonating  : 
the  dynamic  load  becomes  smoother  along  the  contact  course,  in 
value  is  less  than  the  static  load  (if  the  contact  rat : 

The  variation  of  dynamic  load  in  this  region  is  out 
change  of  the  teeth  stiffness  (rig  3-14). 

3.4.2  The  Effect  of  Shaft  Misalignment 

When  the  assembly  errors  are  in. t:  • 


contact  bearing  shifted  to  either 
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Fig.  3-15  shows  the  typical  paths  of  central  contact,  toe  contact,  and 
heel  contact.  Usually  the  central  contact  is  desired  because  it  can 
tolerate  more  possible  running  position  errors  and  avoid  edge  contact. 

The  dynamic  load  response  of  toe  and  heel  contactsare  shown  in  Fig. 

3-16  and  Fig.  3-17.  The  change  of  the  contact  bearing  from  the  center 
to  either  edge  will  also  change  the  contact  ratio  of  the  system  because 
the  tooth  surface  is  not  a  perfect  involute  along  the  profilewise 
direction  and  is  mismatched  along  the  lengthwise  direction.  In  the 
current  example,  the  contact  ratio  for  the  toe  contact  is  1.26, 
central  contact  1.16,  and  heel  contact  1.0.  In  this  case,  if  the  con¬ 
tact  bearing  is  moved  farther  toward  the  heel  region,  there  would  be  no 
tooth  contact  between  the  time  when  the  previous  tooth  finishes  the 
contact  and  current  tooth  goes  into  the  contact  zone  (discontinuity 
in  tooth  mesh).  This  situation  would  cause  a  very  large  impact  force 
which  would  generate  noise  and  severe  damage  to  the  tooth  surface. 

The  effect  of  the  tooth  contact  ratio  on  dynamic  response  is  shown  later. 

3.4.3  Contact  Path  Variation  Due  to  Dynamic  Response 

In  addition  to  showing  the  contact  paths  due  to  assembly  errors 
in  the  system  in  Fig.  3-15,  the  real  contact  path,  including  the  effects 
of  both  assembly  errors  and  running  position  errors  induced  by  the 
dynamic  responses,  is  plotted  in  the  same  figure.  When 
this  real  contact  path  is  compared  with  that  caused  by  the  assembly 
errors  and  running  position  errors  induced  by  the  average  static  elastic 
deformations,  the  deviation  is  found  to  be  surprisingly  small.  One 
explanation  of  this  small  difference  might  be  that  the  displacements 
change  due  to  the  dynamic  oscillation  are  small,  and  they  do  not  produce 
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a  large  change  in  contact  path  compared  to  those  caused  by  the  static 
displacement  only.  The  closeness  between  these  two  contact  paths  suggests 
that  one  can  use  the  average  static  elastic  deformation  to  calculate  the 
contact  path,  which  can  then  be  used  directly  to  solve  for  the  dynamic  load 
and  lubrication  behavior  without  having  to  solve  the  dynamic  load 
and  contact  path  simultaneously  using  an  iterative  technique.  The 
elimination  of  this  iterative  procedure  greatly  reduces  the  computation 
time. 

3.4.4  Effect  of  Speed 

Once  the  physical  conditions  of  a  gear  set  are  determined,  the 
dynamic  response  depends  on  the  operating  speed.  For  a  system  with  one 
degree  of  freedom,  such  as  spur  gears,  the  maximum  dynamic  load  occurs 
when  the  meshing  frequency,  which  depends  on  the  operating  speed,  is 
near  the  system  natural  frequency.  Some  peaks  in  the  dynamic  load 
are  caused  by  the  varied  meshing  stiffness  along  the  contact  path,  and 
they  appear  at  meshing  frequencies  lower  than  the  system  natural  fre¬ 
quency.  The  dynamic  load  factor, defined  as  the  ratio  of  maximum 
dynamic  load  to  the  average  static  load, is  plotted  against  the  gear 
speed  to  illustrate  the  effect  of  speed  in  Fig.  3-18.  Since  there  are 
eleven  degrees  of  freedom  in  the  spiral  bevel  gear  system,  more  peaks 
in  the  dynamic  load  are  expected. 

It  is  seen  that  the  highest  dynamic  load  appears  to  occur  near 
the  natural  frequencies  which  correspond  to  the  mode  associated  with 
the  largest  displacements  in  the  motion  along  the  line  of  action.  The 
frequencies  marked  t  in  Fig.  3-18  shows  the  system  natural  frequency 
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causing  a  larger  displacement  in  the  motion  along  the  line  of  action  and 
the  ones  marked  +  shows  the  system  natural  frequency  with  a  small 
displacement  in  that  motion.  It  is  clearly  shown  that  the  dynamic  load 
factor  at  the  frequencies  marked  *  has  a  peak  response,  and  the  response 
at  the  natural  frequencies  marked  I  is  not  necessarily  a  peak. 

3.4.5  Effect  of  Contact  Ratio 

The  contact  ratio  is  defined  as  the  ratio  of  the  time  required  for 
one  tooth  to  go  through  the  whole  contact  zone  to  the  time  required  for  a 
periodic  meshing  cycle.  It  is  believed  that  the  load  sharing 
characteristics  caused  by  more  than  one  tooth  in  contact  will  reduce  the 
static  load.  The  variations  in  the  dynamic  load  factor  due  to  the  effect 
of  changing  contact  ratio  is  shown  in  Figs.  3-19(a)  to  (c).  It  can  be 
seen  that  the  maximum  dynamic  load  factor  does  not  change  much.  However, 
the  contact  ratio's  effect  is  significant  in  high  speed  regions,  where 
the  load  is  spread  out  between  meshing  teeth  pairs.  A  typical  dynamic 
load  variation  with  a  high  contact  ratio  along  the  contact  path  is  shown 
in  Fig.  3-20. 

3.4.6  Effect  of  Damping 

Since  the  damping  forces  are  usually  not  known  in  the  gear  system, 
three  arbitrary  values  are  chosen  for  the  damping  coefficients:  2627, 
4378,  and  6129  N  sec/m  (15,  25,  and  35  Ib.sec/in).  These  values  are 
selected  to  give  a  range  of  non-dimensional  damping  ratios  corresponding 
to  those  commonly  used  in  spur  gears  (0.1-0. 2).  The  non-dimensional 
damping  ratios  corresponding  to  the  above  damping  coefficients  are  0.087, 


96 


0.14  and  0.203.  The  dynamic  response  for  these  damping  cases  can  be 
observed  from  Figs.  3-18,  3-19(a),  and  3-21.  It  is  expected  that  the 
larger  the  damping  force  is,  the  smaller  the  dynamic  load  factor  will  be 
in  the  resonance  region.  The  large  damping  force  will  also  level  off  the 
peak  of  dynamic  load  factor  in  the  subresonance  region,  and  there  is  no 
effect  on  the  dynamic  load  factor  due  to  the  damping  force  in  the 
superresonance  region. 

3.4.7  Effect  of  Bearing  Stiffness 

It  is  well  known  that  the  bearing  stiffness  plays  an  important 
part  in  the  dynamic  load  response  because  it  directly  controls  the  static 
displacements  which  determine  the  contact  path.  A  large  stiffness  for 
supporting  bearings  is  sometimes  desirable  because  it  pushes  all  the 
resonant  frequencies  beyond  the  range  of  the  operating  speeds.  In  Fig. 
3-22,  the  dynamic  load  factor  is  calculated  for  a  system  with  "infinite" 
bearing  stiffness.  The  contact  ratio  for  this  case  is  1.16,  and  the 
damping  ratios  are  0.058  and  0.14.  Since  the  bearing  stiffnesses  are 
infinite,  only  the  rotational  mode  prevails.  Three  peaks  are  shown  to 
exist  in  Fig.  3-22;  two  of  these  are  the  subharmonics.  These  results  for 
infinite  stiffness  are  compared  with  the  results  shown  in  Fig.  3-19(a) 
for  a  finite  stiffness  of  3.5  *  10**  N/m  (2  x  10^  lb/in.).  It  is  seen 
that  the  natural  frequency  of  the  rotational  mode  with  an  infinite 
stiffness  is  increased  to  290  rad. /sec.  from  the  natural  frequency  of  200 

g 

rad. /sec.  for  a  finite  stiffness  of  3.5  x  10  N/m. 


CONTACT  PATH 


Fig.  3-1.  Dynamic  Modeling  and  Zone  of  Action  of  Bevel  Gears. 


Fig.  5-4.  A  Typical  Section  of  Elements  of  Gear 
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Central  Tooth  Elements  And  Attached  Kira  Elements  of 


Fig.  3-6.  Parts  of  Gea: 


Fig.  3-7.  The  Elements  of  Rim  And  Three  Adjacent  Teeth  of 
Pinion 


beginning  of  contact  end  of  contact 

Fig.  3-9  variation  of  effective  radius  of  pinion 
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Fit,.  3-10.  Stiffness  Variation  Along  The  Contact  path 


Gear  Speed  50  rad/sec  For  Central  Contact 


50  rad.  /sec. 


Speed  150  rad/sec  For  Central  Contact 
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Fig.  3-13.  Dynamic  Load  Variation  Along  The  Contact  Path  at 
Gear  Speed  300  rad/aec  For  Central  Contact. 
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Fig.  3-14.  Dynamic  Load  Variation  Along  The  Contact  Path  at 
Gear  Speed  523  rad/sec  For  Central  Contact. 
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Fig.  3-17.  Dynamic  Load  Variation  Along  The  Contact  Path  at 
Gear  Speed  523  rad/sec  For  Heel  Contact 


Dynamic  Load  Factor  is  Plotted  Against  Gear  Speed 
With  All  Damping  Coefficients  =  15  lbf  •  s/in 
Contact  Ratio  =  1.16 


'Speed  With  Contact  Ratio 


Speed  With  Contact  Ratio 


Fig.  5-20.  Dynamic  Load  Variation  With  Contact  Ratio 


mic  Load  Factoi’  is  Plotted  Against  Gear 
d  With  All  Damping  Coefficients  =  35  lb 
Contact  Ratio  =  1.16 


CHAPTER  IV 


LUBRICATION  OF  SPIRAL  BEVEL  GEARS 


4.1  Introduction 

One  of  the  most  important  failure  modes  in  lubricated  machine 
elements  is  scuffing,  which  is  where  local  wear  occurs  between  the 
contacting  elements  when  the  lubricant  fails  to  separate  them  and  the 
protective  thin  surface  film  breaks  down  due  to  the  high  temperature.  In 
gearing  systems  this  failure  is  affected  mainly  by  gear  geometry,  speed, 
load,  and  lubrication.  For  many  years,  the  methods  used  to  predict  the 
surface  capacity  of  gear  systems  depended  on  empirical  formulas  based  on 
field  experience.  A  major  drawback  in  the  application  of  the  empirical 
methods  is  that  they  do  not  consider  the  effects  of  lubrication,  which 
has  been  found  to  be  of  great  importance  in  scuffing  prediction.  In 
recent  years,  the  fully  developed  theory  in  elastohydrodynamic 
lubrication  was  successfully  applied  to  spur  and  helical  gears  (9,52). 
Because  of  the  geometric  complexity  of  spiral  bevel  gears,  an  analysis 
covering  all  the  possible  factors  to  evaluate  the  risk  of  scuffing  in 
lubricated  spiral  bevel  gear  drives  has  yet  to  be  developed. 

In  this  chapter,  an  analysis  of  the  temperature  and  film  thickness 
for  spiral  bevel  gear  teeth  is  performed.  The  bulk  temperature  is  calcu¬ 
lated  by  using  the  3-D  finite  element  method  combined  with  a  prescribed 
heat  input.  For  the  heavily  loaded  system,  the  limiting  shear  stress 
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depending  on  temperature  and  pressure  is  used  to  calculate  the  total 
heat  generated  over  the  whole  contact  ellipse.  The  average  and  maximum 
flash  temperature  in  this  contact  ellipse  is  obtained  by  making  a 
finite  number  of  stripes.  Each  stripe  is  oriented  along  the  direction 
of  sliding  velocity.  It  is  assumed  that  the  temperature  profile  along 
any  stripe  is  the  same  as  that  of  an  infinitely  long  band  heat  source 
(in  the  direction  perpendicular  to  sliding)  whose  width  is  equal  to  the 
stripe  length  and  has  the  same  heat  flux  profile  along  the  stripe.  The 
heat  flux  partition  coefficient  is  calculated  by  assuming  that  the  total 
heat  is  generated  in  the  mid-plane  of  the  lubricant  film,  and  that  there 
is  no  temperature  jump  in  this  plane.  Hamrock-Dowson' s  equation  is  used 
to  predict  the  film  thickness  between  matching  teeth.  All  quantities 
mentioned  above,  which  include  bulk  temperature,  flash  temperature, 
traction  coefficient  and  film  thickness,  are  interdependent.  They 
are  solved  as  a  system  by  an  iterative  method. 


During  meshing , each  tooth  face  will  show  a  temperature  increase 
(flash  temperature)  due  to  the  frictional  heat  developed  at  the  surface 
as  the  contact  area  moves  along  the  tooth  face.  This  temperature  rise  is 
restricted  to  the  instantaneous  contact  area  and  will  disappear  very 
rapidly  as  soon  as  the  area  of  tooth  face  is  out  of  contact.  Usually, 
this  temperature  is  very  high  and  it  is  an  important  factor  in  gear 
scuffing. 

The  first  successful  prediction  of  flash  temperature,  done 
by  Blok  (13,14),  was  based  on  the  heat  conduction  analysis  of  a 
semi- inf inite  body  with  a  uniformly  distributed  moving  heat  source. 

J.  C.  Jaeger  in  1942  (15)  dealt  with  the  problem  of  a  moving  source  of 
heat  with  variable  heat  source  and  variable  velocity.  J.  F.  Archard 
(16)  introduced  a  simple  harmonic  mean  to  obtain  the  interface  temperature. 
A  refined  solution  including  a  local  heat  partition  coefficient  between 
a  pair  of  disks  was  derived  by  Cameron,  et  al.  (17).  More  recently, 

Francis  (18)  made  a  further  refinement  in  Blok's  calculation  by  con¬ 
sidering  a  variable  heat  flux  in  the  contact. 

The  calculation  of  flash  temperature  is  based  on  the  assumption 
that  the  heat  source  passes  over  the  surface  of  a  semi- inf inite  solid. 

For  spiral  bevel  gears,  the  area  of  sliding  contact  is  formed  by  elastic 
deformation  of  curved  tooth  surfaces  where,  according  to  the  classic 
equation  of  Hertz  (19),  the  area  of  contact  is  elliptical.  Archard  (16) 
has  shown  that,when  the  dimensionless  parameter  vR/a  >  vlO.O  (where  a 
is  the  thermal  diffusivity,  v  is  the  sliding  velocity  and  R  is  the 
nominal  length),  the  heat  flow  in  the  direction  perpendicular  to  sliding 
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may  be  neglected.  The  temperature  distribution  within  a  heat  source 
of  finite  area  can  be  determined  by  dividing  the  whole  contact  area 
into  differential  stripes  parallel  to  the  sliding  direction.  The 
temperature  profile  along  any  stripe  is  the  same  as  that  of  an  in¬ 
finitely  long  band  heat  source  (in  the  direction  perpendicular  to  slid¬ 
ing)  of  width  equal  to  the  stripe  length  and  having  the  same  heat  flux 
profile  along  the  stripe. 

The  situation  is  illustrated  in  Fig.  4-1,  which  shows  the 

contact  heat  flux  distribution  in  an  elliptical  contact  area  with  an 

angle  a  between  direction  of  the  velocity  of  heat  flux  (sliding  velocity) 

2 

and  semi-minor  axis  of  this  area.  Letting  the  equation  (X/AMAX)  + 

2 

(Y/BMAX)  =  1  describe  this  elliptic  contact,  the  length  g  of  a  stripe  con¬ 
taining  point  (x,y)  in  the  contact  is 

YA  =  BMAX2  tan2  0  +  AMAX 
YB  =  tan  0  BMAX2  x^  -  y  BMAX2  tan2  0 

YC  =  BMAX2  tan2  0  y2  -  2  tan  0  BMAX2  x  y± 

+  BMAX2  x2  -  AMAX2  BMAX2 

+  ,-  -  YB  ±  VyB2  -  YA  YC 

y  a  - 

YA 

*+•"  =  tan  0  (y+’~  ~  yi)  +  xi  (4.5) 


(4.1) 

(4.2) 

(4.3) 

(4.4) 
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21  =  ^(y*  -  y")2  +  (X*  -  X-)2  (4.6) 

The  center  of  this  stripe  is 

xiQ  =  0.5  (x+  +  x±)  (4.7) 


y.  =  0.5 
10 


(yi  +  V 


(4.8) 


The  dimensionless  parameter  u  is  defined  as 


u 


x. 

lO 


)2  +(y  -  ylo> 


l 


(4.9) 


If  (x,y)  is  in  the  forward  position  of  this  stripe  (upstream  of  the 
sliding  velocity),  u  is  positive;  otherwise,  u  is  negative.  The 
temperature  at  this  point  can  be  obtained  by  using  the  equation  developed 
by  Jaeger  (15)  for  the  temperature  distribution  along  a  stripe  for  a 
farger  value  of  v£/a. 


T(x,y)  = 


-  IS-  I  a(ft-u) 


rv 


1/2 


(4.10) 


where  k  is  the  thermal  conductivity  and  q  is  a  constant  heat  flux 
(the  reason  for  using  a  constant  heat  flux  will  be  shown  later)  Cameron, 
et  al.  (17)  showed  that,  for  a  uniform  source,  the  temperature  profile 
is  given  almost  exactly  for  10  <  v£/a  <  ®  by  multiplying  the  asymptotic 
(vi/a  'v*  00 )  factor  1  +  (2i/i-u)^ ’ ^ [0. 65  (a/vi)  +  0.44  (a/vi)^^2],  and 
the  temperature  profile  becomes 
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T(x,y)  =  & 


a(£-u)l^'3 

f  Oj?  'I  ®  ^  3/2 

1  +  hrH  [0. 65 (a/v£)  +  0.44(a/vO  7  ]' 

irv  J  I 

l  Jc— ul 

i 

(4.11) 


This  factor  will  increase  the  maximum  temperature  and  moves  it  toward  the 
center  of  the  stripes  as  vl/a  decreases.  Even  if  vi/a  <  10,  which  occurs 
in  the  outer  stripes  of  the  contact  ellipse,  the  equation  can  still  be 
used  because  (1)  the  correction  factor  will  hold  approximately  for 
vl/a  <  10.0,  and  (2)  the  outer  stripes  are  only  a  small  fraction  of  the 
area  of  the  heat  source.  Thus,  there  will  be  a  minimal  effect  on  the 
temperature  in  the  whole  contact  area. 

Heat  is  generated  by  viscous  shearing  of  the  lubricant  between  the 
sliding  surfaces.  This  heat  is  carried  away  either  by  the  lubricant 
through  convection  or  by  the  gear  blanks  through  conduction.  The 
relative  importance  between  these  two  modes  of  heat  transfer  in  EHD 
contact  was  examined  by  Trachman  (20).  He  concluded  that  the  heat 
carried  away  by  heat  convection  can  be  neglected  in  comparison  to  the 
heat  diffused  to  the  gear  body,  except  at  very  high  speeds  which  exceed 
those  of  gears  in  current  practice.  Let  the  subscripts  1  and  2  denote 
the  pinion  and  the  gear,  respectively.  The  heat  generated  in  the  contact 
area  is  divided  betwen  the  pinion  and  the  gear  by  the  ratio  a/(l  -  a), 
where  a  is  the  heat  partition  coefficient.  Therefore,  within  the 
contact,  the  temperature  distribution  in  the  pinion  and  the  gear  are 


T  (x,y)  =  a 


(WV 

0.5 

1  + 

'  Zl  ' 

TTV- 

< 

X  ' 

£  ,-u. 

1 

1  1 

\ 

1 

0.5 


[0.65  (a^/v^) 


+  0.44  a1/v1^1)3/2]l 


(4.12) 
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T2(x,y) 


(1-a) 


23. 


a2(£2-u2) 


irv„ 


\  0.5 

! 

2. 

< 

1  + 

2 

lV“2. 

0.5 


[0.65(a2/v2£2) 


+  0.44(a2/v2£2)3/2]' 


(4.13) 


The  lubricant  local  velocity  profile  and  temperature  profile 
for  the  very  high  sliding  velocity  case  were  estimated  by  Plint  (21) 
and  later  confirmed  by  Trachman's  analysis  (20).  For  sliding  contacts, 
their  results  showed  a  sharp  S-shape  velocity  profile  across  the  film 
with  a  large  velocity  gradient  at  the  mid-plane.  The  temperature  is  at 
its  maximum  at  the  mid-plane  and  decreases  almost  linearly  to  both 
surfaces.  Since  shearing  of  the  lubricant  occurs  mainly  in  the  mid¬ 
plane,  most  of  the  heat  is  also  generated  in  this  layer.  This 
results  in  a  triangular  temperature  distribution  across  the  film. 

According  to  Francis'  analysis  (18),  if  the  bulk  surface 
temperatures  for  the  pinion  and  the  gear,  TR^  and  are  different 

(see  Fig.  4-2),  then  the  interfacial  temperature  can  be  expressed 


(a  +  Aa) 

r  z  m 

ti+  ; 

O 

v. 

I 


-  (1-a-Aa) 


+  T. 


B1 


mid 


(x,y) 


t2  + 


(h-Zjq 

m 


+  T 


B2 


(4.14) 


where  Aa  is  the  unknown  function  of  (x,y)  which  expresses  the  local 

deviation  from  heat  flux  distribution  aq  and  (1  -  a)q,  Z  is  the 

m 

distance  from  pinion  surface  to  the  layers  at  which  all  the  shear 


is  concentrated  (Z  =  — ) ,  k  is  the  oil  thermal  conductivity,  and 

m  z  o 

T  ^(x.y)  is  the  temperature  at  the  shear  plane.  Then 
mid 


(a  +  Aa) 
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T  + 

1  2k, 


+  TBi  =  C1  “  «  “  Aa) 


hq 

T  +  — — 
2  2k 


+  T. 
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(4.15) 


A  +  Act  = 


T2  +  hq/2ko  +  Tb2  -  Tb1 

T,  +  T,  +  hq/k 
1  ^  o 


(4.16) 


Tmid(x’y) 


T2(x,y)  +  hg/2kQ  4-  Tb2  -  TB1 

T  (x,y)  +  T,(x,y)  +  hq/k 
l  c  o 


( 

Tx(x,y) 


+  T 


Bl 


(4.17) 


The  heat  partition  coefficient  a  must  now  be  determined  to  calcu¬ 
late  the  flash  temperature  of  pinion  and  gear.  The  analytical  expres¬ 
sion  for  a  in  terms  of  local  conditionals  difficult.  The  best  approxima¬ 
tion  of  a  will  be  the  value  that  satisfies  the  condition  for  which  the 
temperature  fields  of  the  two  surfaces  are  identical  over  the  contact 
ellipse.  Thus,  an  analytical  approximation  can  be  written  by  equating 
the  mean  temperature  of  the  two  surfaces,  as  done  by  Jaeger  (15) 


r  1 

r  i 

Qh 

T  +  — — 

+  =  (1  -  a) 

Qh 

To  + 

l,ave  2k 

o; 

Bl 

2,ave  2k 

o 

+  T 


B2 


(4.18) 


where 
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If  T1(x,y)dxdy 


l.ave  it  AMAX  BMAX 


(4.19) 


If  T2(x,y)dxdy 


2,ave  tt  AMAX  BMAX 


(4.20) 


Q  is  the  total  heat  flux  over  the  whole  contact  ellipse.  Once 
the  a  is  determined,  the  flash  temperatures  become 


Il  f<x,y)  -  aTjh.y)  +  TB1 


(4.21) 


T2  f(x,y)  =  (1-a)  T2(x,y)  +  Tb2 


(4.22) 


129 


4.3  Bulk  Temperature 

The  method  mentioned  in  Section  4.2  is  for  calculating,  with  known 
heat  flux,  the  flash  temperature  rise  above  the  surface  temperature  which 
was  assumed  known  beforehand.  The  existing  criteria  for  scuffing  of 
heavily  loaded  gears  is  that  the  total  surface  temperature  rises  to  a 
critical  value  at  which  the  oil  film  collapses,  the  protective  surface 
layer  breaks  and  metal  to  metal  contact  occurs.  Uncertainty  in  the  bulk 
surface  temperature  makes  any  scuffing  prediction  useless  even  if  an 
accurate  calculation  of  flash  temperature  can  be  calculated. 

Before  the  gear  system  starts  to  operate,  all  elements  are  at  the 
ambient  temperature.  Once  the  gear  system  starts  to  operate,  the 
temperature  builds  up  as  the  gears  are  running  due  to  frictional  heat. 
After  a  sufficient  period,  the  gear  surface  temperatures  reach  a  steady 
state  at  which  the  heat  flux  flowing  into  the  body  is  equal  to  that 
flowing  out  of  the  body.  At  each  revolution,  the  tooth  is  subject  to  the 
same  heating  flux.  Since  the  time  period  for  each  contact  point  in  the 
contact  zone  is  only  a  small  fraction  of  the  entire  period  of  revolution, 
the  local  temperature  jump  (flash  temperature)  has  completely  decayed 
before  it  enters  into  the  contact  zone  at  the  next  revolution.  Thus,  it 
is  justified  to  use  an  average  heat  input  over  the  revolution  to 
calculate  the  steadjrstate  temperature  rise  of  the  body. 

The  heat  input  is  due  to  the  heat  generated  at  the  instantaneous 
contact  ellipse,  and  the  amount  depends  on  the  load  and  the  viscous 
shearing  of  the  lubricant.  The  heat  flux  flowing  out  of  the  body  is  due 
to  the  heat  convection  to  the  surrounding  air  and  lubricant.  The 
relative  importance  of  the  heat  transfer  coefficient  at  different  surface 
areas  was  discussed  by  Patir  (22)  and  Townsend  and  Akin  (23)  in  spur  gear 
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systems.  They  also  revealed  the  significant  effect  on  temperature 
distribution  by  the  lubrication  method.  In  this  study,  the  oil  jet 
impingement  depth  is  assumed  to  cover  the  whole  area  of  contact  side, 
which  can  be  obtained  by  using  the  proper  pressurized  oil  jet.  Estimated 
values  of  the  heat  transfer  coefficients  are  used  on  the  other  surface 
areas  to  calculate  the  bulk  temperature. 

A  three-dimensional  finite  element  program  is  used  to  calculate 
the  temperature  coefficient.  The  mesh  of  the  system  includes  the  gear 
shaft,  gear  body  and  the  contact  tooth  with  one  adjacent  tooth  on  both 
sides.  The  eight-node  element  is  used.  This  model  is  the  same  as  that 
used  for  the  deflection  coefficient  except  that  the  boundary  conditions 
are  different.  In  the  temperature  analysis,  all  the  surfaces  are  subject 
to  heat  convection  with  different  heat  transfer  coefficients  except  for 
the  surface  in  the  inner  cylinder  of  the  shaft,  which  is  assumed  to  have 
no  heat  convection.  The  inner  cylinder  of  the  gear,  created  for  the 
modeling,  is  very  small  in  radius,  therefore  heat  transfer  can  be 


neglected  for  these  surfaces.  Since  the  air  in  the  inner  cylinder  of  the 
pinion  shaft  is  enclosed,  there  is  no  heat  loss  or  gain  in  the  steady 
state,  and  the  heat  convection  at  this  surface  can  also  be  neglected. 

The  heat  transfer  coefficient  hj  is  assigned  to  the  contacting  tooth  face 
which  is  oil  jet  cooled.  The  topland  and  bottomland  of  the  tooth  and  the 
other  side  of  the  tooth  surface,  which  are  not  cooled  by  the  oil  jet, 
will  have  a  heat  transfer  coefficient  ht  for  air  or  air/oil  mist.  Since 
only  three  teeth  are  made  in  the  model,  there  is  a  surface  region  A 
covering  the  surfaces  where  the  teeth  are  taken  off  and  the  bottomland 
between  those  teeth.  The  heat  transfer  coefficient  at  this  region  A  is 


given  a  value  of  h^  ,  which  is  the  same  as  that  of  the  surface  coefficient 

cooled  by  the  oil  jet.  The  reason  is  that  there  is  an  oil  jet  cooled 

surface  on  each  tooth  and  most  of  the  heat  will  flow  out  of  the  tooth  through 

this  surface  (h^  >>  ht).  All  the  other  convective  surfaces  of  the  gear 

system  are  given  a  coefficient  hg  (Fig.  4.3).  The  theoretical  estimated 

values  of  h  and  h.  can  be  found  in  (Refs.  24  and  25).  However,  the 
s  j 

estimated  values  of  hg ,  h  and  h^  based  on  the  experimental  results  are 
used  in  this  study  (23). 

Thirty  nodes  are  created  on  the  contacting  surface.  The  tempera¬ 
ture  distribution  on  this  surface  is  represented  by  T  which  is  the 
temperature  at  the  grid  node  J  due  to  a  unit  heat  flux  at  node  I.  By 
interpolation,  the  temperature  at  any  contact  point  M  due  to  a  unit  heat 

Q 

flux  at  the  contact  point  N  (T  T)  can  be  obtained  in  terms  of  Ttt.  The 

MN  J1 

details  of  this  interpolation  can  be  found  in  Appendix  D.  Once  the 
contact  path  is  located  and  the  heat  flux  flowing  into  each  body  at  each 
contact  point  has  been  calculated,  the  bulk  temperature  at  the  contact 
point  M  can  be  found  as 

TM,B=ZtCMNQN  N=1...KMAX  (4.23) 

where  KMAX  is  the  total  number  of  contact  points  along  the  contact  path. 
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4.4  Traccion  Coefficient 


In  analyzing  lubrication  related  failure  mode  in  EHD  contacts,  the 


heat  generated  by  traction  in  the  contact  area  is  the  main  concern.  In  a 


conventional  analysis,  the  rheology  of  the  lubricant  is  ignored  and  the 


local  shear  stress  in  the  contact  area  is  assumed  to  be  proportional  to 


the  local  pressure  and  acts  in  the  direction  opposite  to  the  local  slip. 


However,  recent  research  into  the  shear  behavior  of  EHD  makes  it  possible 


to  examine  the  role  of  the  fluid  rheology. 


The  magnitude  of  the  coefficient  of  friction  in  a  highly  loaded 


disk  machine  was  examined  by  Misharin  (26)  at  high  rolling  and  sliding 


speeds.  He  also  gave  an  approximated  formula  for  the  prediction  of  the 


friction  coefficient  based  on  his  experimental  data.  But  his  formula  did 


not  give  realistic  results  at  a  very  high  or  very  low  value  of  the 


rolling  or  sliding  speeds.  Smith  (27)  presented  a  set  of  results  that  showed  that 


the  friction  coefficient  would  increase  to  a  maximum  and  then  decrease 


as  the  sliding  speed  is  increased.  He  pointed  out  that  the  decrease  of  the 


friction  coefficient  was  due  to  thermal  and  non-Newtonian  effects. 


The  thermal  effect  upon  the  viscosity  of  lubricant  was  examined 


by  Crook  (28,29)  and  by  Cheng  (30).  Crook  used  a  disk  machine  to  study 


rolling  friction  and  sliding  friction.  He  observed  that  the  plot  of 


traction  against  slip  always  showed  a  same  kind  of  characteristic  pattern 


as  that  which  was  found  by  Smith  (27)  regardless  of  the  variation  of  load 


and  speed.  He  also  attempted  to  predict  the  traction  analytically  by 


simplified  friction  theory  based  on  the  following  assumptions:  The  film 


thickness  in  the  contact  area  is  uniform;  the  pressure  distribution  in  the 


contact  area  is  Hertzian;  the  heat  carried  away  by  convection  through  the  lubricant 
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is  neglected;  and  the  temperature  rise  on  the  surface  is  neglected.  His 


method  can  predict  a  traction  coefficient  consistent  with  experimental 


data  in  high  sliding  region.  However,  he  could  not  explain  the  behavior 


in  low  sliding  region.  Cheng  used  a  thermal  EHD  theory  which  included 


energy  equation  and  was  free  of  all  the  assumptions  made  by  Crook. 


Cheng's  results,  compared  with  Crook's  experimental  data,  still  showed 


a  big  difference  in  friction  coefficient  in  the  low  sliding  region.  Both 


works  suggested  that  the  thermal  effect  alone  cannot  account  for  the 


sharp  reduction  of  the  effective  viscosity  in  the  low  sliding  region. 


A  pioneering  effort  was  made  by  Dyson  (31).  He  studied  various 


experimental  data  and  found  that  the  curve  of  traction  coefficient,  when 


plotted  against  sliding  speed,  can  be  divided  into  three  regions  (Fig. 


4-4).  Region  I  is  the  linear  region  which  shows  the  characteristic 


behavior  of  a  Newtonian  fluid.  Region  II  is  the  nonlinear  region  which 


shows  nonlinear  variation  of  the  traction  coefficient  with  the  shear  rate,  and 


region  III  is  the  thermal  region  which  is  dominated  by  the  thermal  effect 


on  the  shear  stress  of  the  lubricant.  The  fluid  behavior  in  region  I 


is  shown  to  be  explainable  on  the  basis  of  linear  viscoelasticity.  In 


region  II,  the  viscosity  will  decrease  when  the  material  undergoes  a 


steady  continuous  shear  with  large  strain.  Dyson  (31)  and  Gruber  and 


Litvitz  (32)  suggested  that  the  nonlinear  behavior  of  the  traction  curve 


is  due  to  this  effect.  At  a  high  pressure,  the  relationship  between 


stress  and  strain  for  the  fluid  is  nonlinear;  several  equations  were 


developed  to  describe  this  nonlinear  behavior.  Bell,  Kannel  and  Allen 


(33),  Hirst  and  Moore  (34)  used  a  "sinh"  relation,  and  Trachman  and  Cheng 


(35)  used  a  hyperbolic  model.  Observing  that  the  traction  force  never 
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seems  to  exceed  one  tenth  of  the  normal  load.  Smith  (36),  Johnson  and 
Cameron  (37),  and  Plint  (38)  proposed  that  the  lubricant  film  would  fail 
to  retain  a  very  high  shear  stress  at  extremely  large  strain  rate,  and 
shears  like  a  plastic  solid. 

More  recent  work  on  the  rheological  experiment  was  done  by 
Johnson  and  Roberts  (39).  They  observed  that  the  liquid-solid  transition 
was  occurring  in  EHD  contact  and  the  transition  point  was  dependent  on 
temperature  and  pressure.  The  same  observations  were  made  by  Johnson 
and  Tevaarwerk  (40)  and  Bair  and  Winer  (41,42).  Johnson  and  Tevaarwerk 
(43)  then  developed  an  elastic-perfect  plastic  solid  model  for  the  traction 
drive  analysis.  Based  on  their  experimental  results,  most  recently 
Bair  and  Winner  (44)  proposed  the  following  model  to  predict  the  lubri¬ 
cant  behavior  over  the  whole  range  of  visco-elastic-plastic  change. 

4.4.1  Bair-Winner  Model 

Bair  and  Winner  fabricated  an  apparatus  to  measure  the  viscosity  of 
thelubricant  at  high  strain  rates  and  a  large  shear  stress  similar  to 
those  in  an  EHD  contact.  They  found  that  the  lubricant  under  high  pres¬ 
sure  exhibits  a  classical  visco-elastic  behavior  for  small  strain.  For 
large  strain  and  large  stress,  the  lubricant  behaves  like  a  material 
having  a  limiting  yield  shear  stress.  For  lubricant  under  a  more 
moderate  pressure  but  large  strain,  the  result  is  that  the  shear  stress  is 
smaller  and  the  lubricant  is  viscous. 

From  their  experimental  results  for  5P4E  Santotrac-50  at  large 
strain,  the  relationship  of  the  dimensionless  shear  stress  and  the 
dimensionless  shear  rate  can  be  described  reasonably  well  with  a  single 
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natural  log  function.  The  dimensionless  shear  stress  t  is  defined  as 


the  actual  shear  stress  t(P,T)  divided  by  the  limiting  shear  stress 


tl(P,T),  and  the  dimensionless  shear  rate  y  *s  the  actual  shear  rate  y 


multiplying  by  the  viscosity  y  (P,T)  and  divided  by  the  limiting  shear 

o 


stress  x  (P,T).  t(P,T),  t  (P,T)  and  y  (P,T)  are  all  pressure  and 

L  L  O 


temperature  dependent. 


This  natural  log  function  can  be  expressed  by 


Y  ■  -  £n(l  -  x) 


(4.24) 


which  governs  the  nonlinear  viscous  flow. 


Based  on  the  conventional  Maxwell  visco-elastic  model,  the  total 


shear  rate  of  the  lubricant  consists  of  the  shear  rate  of  elastic  part 


(Ye)  and  that  of  the  viscous  part  (y^) 


y  -  y  +  y 
'  'e  v 


(4.25 


The  elastic  part  can  be  nondimensionalized  the  same  way  as 


Y  M  (P.T)  x(P,T)  y  (P.T) 


tl(P,T) 


G  (P,T)t  (P.T) 
00  L 


where  G  (P,T)  is  the  limiting  elastic  shear  modulus  at  a  given  tempera¬ 


ture  and  pressure.  With  the  above  natural  log  relation  for  the  viscous 


part,  a  modified  Maxwell  model  is  obtained: 


'tVV-.V  V  V  v.v 


r* i> 1 


In  order  to  use  this  equation,  three  primary  physical  properties  are 


required.  These  include  the  low  shear  stress  viscosity  sjq,  the 
limiting  yield  shear  stress  t  ,  and  the  limiting  elastic  shear  modulus 
G^.  These  three  parameters  are  all  functions  of  the  temperature  and 
pressure . 

In  the  application  of  the  spiral  bevel  gears  in  heavy  loaded 
cases,  the  pressure  in  most  contact  ellipses  is  very  high  (1  G  Pa  to  2 
G  Pa  or  150,000  psi  to  300,000  psi)  and  the  slide  to  roll  ratio  is  about 
0.05  -  0.3.  Under  these  severe  conditions,  the  limiting  shear  stress 
for  most  lubricants  should  occur  in  most  of  the  contact  area.  For  the 
temperature  calculation,  it  is  safe  to  assume  that  the  uniform  limiting 
shear  stress  is  applied  at  a  given  average  pressure  and  temperature  in 
the  contact  ellipse. 
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(4.28) 


where  is  given  by  experimental  data. 

Since  experimental  data  for  tt  in  the  temperatures  encountered  in  spiral 

Li 

gears  are  not  available,  the  conventional  way  to  calculate  this  yield  shear 
stress  postulated  by  Dyson  is  used: 

tl  =  GJ  4  (4.29) 

where  G^  is  called  limiting  shear  modulus,  and  G^  is  a  function  of 
temperature  and  pressure.  Based  on  the  experimental  work  of  Hutton  (11) 
and  Switch  (27),  Dyson  proposed  a  relationship  between  G^,  temperature 
and  pressure  for  high  viscosity  index  mineral  oils: 
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A  subroutine  was  prepared  to  evaluate  the  traction  based  on  Bair 
and  Winer's  model.  Unfortunately,  for  the  temperatures  encountered  in 
spiral  bevel  gears,  the  values  of  pQ,  tl  and  are  not  available  from 
Bair  and  Winer's  work. 
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4 .5 _ Film  Thickness 

In  1916,  Martin  provided  the  first  paper  about  roller  lubrication. 
Although  he  failed  to  predict  a  correct  film  thickness  by  using  the 
assumptions  of  a  rigid  body  and  an  incompressible  isoviscous  fluid,  his 
approach  was  welcomed  by  later  investigators.  The  influence  on  the  film 
thickness  by  the  pressure  dependence  of  the  viscosity  in  the  contact  area 
was  examined  by  Gatcombe  (45)  by  using  an  exponential  relationship  and  by 
Hersey  and  Lowdenslager  (46)  by  using  parabolic  viscosity-pressure 
relationships.  More  investigations  in  this  field  were  made  by  Cameron 
(49)  and  by  McEwen  (48).  These  results  showed  that  the  effect  of 
pressure  dependence  on  viscosity  could  not  alone  account  for  the 


difference  between  prediced  and  observed  film  thickness,  although  it  did 
provide  some  improvement.  The  effects  of  elastic  deformation  and 
pressure  dependent  viscosity  were  considered  first  by  Ertel  and  Grubin  in 
1945,  who  initiated  the  study  of  EHD  lubrication.  By  assuming  the 
deformation  in  the  inlet  region  of  contact  zone  was  the  same  as  that  of 
the  dry  Hertzian  contact,  they  calculated  the  pressure  and  film  thickness 
in  this  region.  Dowson  and  Higginson  (19)  developed  a  general  formula 
for  minimum  film  thickness  for  line  contacts  which  included  the  effects 
of  elastic  distortion  and  pressure  dependent  viscosity.  The  dimension¬ 
less  form  of  this  equation  is 


where 


H 


min,  £, 


1.6  G°’6  U0'7  WD°*13 


G  =  oE 

U  =■  n  u/e'r 
o 

wQ  -  w/e'r 


(4.31) 
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R  is  the  effective  radius  of  the  roller  pair,  E'  is  the  effective  elastic 
modulus,  nQ  is  the  atmosphere  viscosity  of  the  lubricant,  a  is  the 
pressure  dependent  viscosity  parameter,  n  =  h  exp(ctp),  w  is  the  load  per 
unit  length,  and  u  is  the  speed.  The  thermal  effect  due  to  sliding 
was  examined  by  Cheng  and  Stemlicht  (50)  by  using  numerical  technique 
to  solve  the  Reynold,  elastic,  energy  and  heat  transfer  equations.  Their 
results  indicated  that  the  isothermal  Dowson-Higginson  equation  could 
still  predict  a  good  approximation  to  the  measured  film  thickness  if  nQ 
was  taken  as  the  value  at  the  bulk  surface  temperature  of  the  contact 
bodies . 

The  'cylinder  model'  used  in  the  lubrication  analysis  of  gears 
is  only  good  at  the  pitch  point.  Actually,  the  effective  radius  and  rolling 
velocity  are  varied  from  point  to  point  along  the  contact  path.  Wayne 
and  Rodzimovsky  (51)  were  among  the  first  investigators  who  examined  the 
actual  contact  ratio,  actual  involute  profile,  the  combined  effects  of 
rolling  and  sliding,  and  film  thickness  in  the  whole  course  of  the 
contact  path.  Unfortunately,  their  analysis  did  not  include  the  elastic 
deformation  in  the  contact  zone.  Based  on  the  quasi-steady  state  assumption, 
Gu  (52)  used  the  EHD  lubrication  theory  to  calculate  the  film  thickness 
of  the  spur  gears  by  including  the  time  dependent  term  3(ph)/3t  in  the 
Reynolds  equation.  He  found  that  this  term  did  not  have  much  effect  on 
the  minimum  film  thickness. 

In  spur  gears,  the  contact  between  mating  teeth  is  a  line 
contact  extended  from  one  end  of  the  tooth  to  the  other  end  of  the  tooth. 
Provided  that  the  bulk  surface  temperature  of  the  gear  teeth  is  known 
beforehand,  one  cannot  predict  accurately  the  minimum  film  thickness  by 


using  the  Dowson-Higginson  equation.  For  conjugate  spiral  bevel  gears,  the 
contact  is  still  a  full  line  from  one  limit  to  the  other  limit.  However, 
fully  conjugated  spiral  bevel  gears  seldom  exist  in  practice.  The  actual 
spiral  bevel  gears  are  crowned  in  both  directions  (lengthwise  and  profile- 
wise)  to  reduce  sensitivity  to  the  shift  of  contact  patterns  due  to 
misalignment.  In  such  cases,  an  elliptical  contact  area  appears  between 
the  mating  teeth  under  load,  for  which  line  contact  solution  is  no  longer 
valid . 


Hamrock  and  Dowson  (53)  evaluated  numerically  the  analysis  of  an 
isothermal  EHD  lubrication  point  contact.  They  showed  that  the  minimum 
film  thickness  could  be  related  to  the  well-known  line  contact  solution 
by  simple  expression  involving  the  ellipticity  parameter  (54).  In  this 
study,  Hamrock-Dowson' s  point  contact  solution  will  be  adapted  in  the 
elliptical  contact  of  spiral  bevel  gears. 

The  effective  situation  of  contact  between  the  spiral  bevel  gear  and 
pinion  can  be  seen  in  Fig.  4-5,  in  which  there  is  a  flat  plane  contact 
with  a  body  which  is  described  by  the  difference  between  neighboring 
surfaces  of  the  gear  and  the  pinion  at  the  contact  point.  This  curved  body 
has  effective  radii  Rx  and  along  the  principal  axis  x  and  y,  respec¬ 
tively.  Under  a  load  P,  the  surface  near  the  flat  plane  will  deform 


to  an  elliptical  shape  with  semi-major  axis  AMAX  and  semi-minor  axis 


BMAX.  V  and  V  are  the  velocities  of  the  pinion  and  gear  at  the  contact 
P  g 

point.  The  ellipticity  parameter  was  defined  as  the  ratio  of  AMAX  to 


BMAX.  The  minimum  film  thickness  in  the  contact  zone  was  related  to 
Dowson-Higginson' s  line  contact  solution  by  the  equation 
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H  =  H  .(1.0  -  1.6  e“°‘62k) 
min  min , l  ' 


(4.32) 


where  H 


^  :  Dimensionless  film  thickness  of  Dowson-Higginson  solution 


K  :  Ellipticity  parameter 
h 


and  H  .  .  “min, l 

min,£  =  — jp1 
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U  =  V  cos  9 
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U  =  V  cos  0 
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U  =  (U  +  U  )/2.0 
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U  -  (U  +  U  )/2.0 

y  py  gy 


=  Vl 


V  =  \  Ux2  +  Uy2 


The  dimensional  form  h  .  .of  the  minimum  film  thickness  for  a  line  contact  can 

min,  & 

be  written  as 


.  .  ,  0.6,  ,n0.7  ,0.43  .  0.13 

h  .  =  1.6  a  (n  V)  E  /w 

min,£  o 


(4.33) 


where  w  *  P/2AMAX  is  the  load  per  unit  length  and  P  is  total  load.  The 
viscosity  nQ  is  taken  at  the  bulk  surface  temperature  to  account  for  the 
thermal  effects  on  the  film  thickness. 


It  is  important  to  note  that  nQ»  P,  AMAX,  BMAX  and  V  are  varied 
along  the  path  of  contact.  V  depends  on  the  gear  kinematics,  and  AMAX 
and  BMAX  depend  on  the  gear  geometry  and  dynamic  load  P.  nQ  is  strongly 
dependent  upon  the  local  static  surface  temperature  which,  in  turn,  is 
influenced  by  the  local  film  thickness  through  frictional  heating.  Thus, 
the  film  thickness  and  the  static  (bulk)  surface  temperature  are 
interdependent,  and  must  be  solved  as  a  coupled  system.  The  dynamic  load 
is  assumed  to  be  uninfluenced  by  the  film  thickness. 


r-  v*  v* *r< 
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4.6 _ Results  of  Lubrication  Performance 

The  same  set  of  gears  used  for  the  dynamic  load  calculations  are 
used  here  to  demonstrate  calculation  of  the  lubrication  performance.  Results 
were  obtained  for  a  range  of  operating  conditions  to  determine  the  effect 
of  speed,  load,  lubricant  viscosity,  and  ambient  temperature  on  the  film 
thickness,  bulk  temperature,  and  flash  temperature. 

The  sliding  velocity  decreases  from  the  beginning  of  the  contact 
path  where  the  gear  tip  contacts  the  pinion  root,  until  the  contact  point 
is  near  the  pitch  point  where  the  sliding  velocity  becomes  zero.  Then  the 
sliding  velocity  increases  all  the  way  to  the  end  of  the  contact  path 
where  the  pinion  tip  contacts  the  gear  root.  The  current  set  of  gears 
has  the  feature  that  the  sliding  velocity  at  the  end  of  the  contact  path 
is  larger  than  that  at  the  beginning  of  the  contact  path;  this  fact 
creates  a  situation  where  more  heat  is  generated  at  the  end  of  the  contact 
path. 

Typical  distributions  of  the  bulk  temperature  along  the  contact 
path  are  shown  in  Figs.  4-6  through  Fig.  4-9  for  various  speeds. 

The  bulk  temperature  of  the  pinion  is  always  larger  than  that  of  the  gear 
because  the  pinion  speed  is  three  times  faster  than  the  gear  speed,  and 
receives  more  heat  in  a  unit  time  than  the  gear  does.  Although  the 
temperature  coefficients  are  higher  near  the  gear  tip,  the  maximum  bulk 
temperature  of  both  gears  occurs  at  the  end  of  the  single  tooth  contact 
point  where  the  maximum  heat  is  generated.  Distributions  of  the  total 
flash  temperature  for  the  same  cases  are  plotted  on  Fig.  4-10  through 
Fig.  4-13.  The  minimum  flash  temperature  occurs  at  the  pitch  point  where 
the  sliding  velocity  is  zero.  For  the  high  speed  case,  the  variation 


of  dynamic  load  is  less  pronounced  along  the  contact  path.  The  rise  of 
flash  temperature  on  both  sides  of  the  pitch  corresponds  directly  to  the 
variation  of  sliding  speed  at  the  contact.  The  slight  decrease  at  the 
end  of  the  contact  path  is  attributable  to  the  decrease  in  dynamic  load 
in  this  region. 

Fig.  4-14  shows  the  distributions  of  the  film  thickness  for  four 
different  gear  speeds.  No  excessive  variations  are  seen  along  the 
contact  path.  A  moderate  peak  is  evident  at  the  pitch  point  for  the  high 

speed  cases,  and  this  is  associated  with  the  slight  drop  of  bulk  tempera¬ 
ture  at  the  pitch  point.  The  steady  rise  of  film  thickness  along  the 
contact  path  is  due  to  the  increase  in  the  entrainment  velocity.  The 
final  uptrend  of  film  thickness  near  the  end  of  contact  is  again  due  to 
the  decrease  in  the  bulk  temperature. 

Finally,  the  effects  of  an  increase  in  ambient  viscosity,  ambient 

temperature,  and  load  on  the  minimum  film  thickness  h  .  ,  maximum  bulk 

min 

and  total  flash  temperatures  is  demonstrated  in  Fig.  4-15.  As 
expected,  an  increase  in  viscosity  would  improve  lubrication  performance 
with  a  much  greater  film  thickness  and  a  slight  drop  in  both  bulk  and 
flash  temperature.  An  increase  in  ambient  temperature  would  reduce 
the  film  thickness  considerably,  and  an  increase  in  load  would  also 
reduce  the  film  thickness  due  to  the  increase  in  the  bulk  temperature. 


body  2 


body  1 


Fig.  4 -2  temperature  distribution  across  film 
thickness  (h)  with  two  different  bulk 
temperatures 


149 


I  \  ' 


& 


"Hi 


<■  i  l 

Vi 


■Hlj-ji 


0 


Vi 


b'fl 


Effective 

radius 


Effective 

radius 


\  / 
\  / 


YgA 


Contact  Situation  And  Contact  Area  With 
of  Pinion  And  Gear 


r-  4%  1 


WWW! 


153 


* 


<v 


« *| 

t  •] 
.  *1 


CONTACT  CONTACT 


Bulk  Temperature  Distribution  For  Gear  Speed 
523  rad/sec 


154 


Fig.  4--10  flash  temperature  distribution  for  gear  speed  50 
rad/sec 


beginning  of  contact  end  0f  contact 


Fig.  4-12  flash  temperature  distribution  for  gear  speed  300 
rad/sec 
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Fig.  4-13.  Flash  Temperature  Distribution  For  Gear  Speed 
523  rad/sec 
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Fig.  4-14.  Minimum  Film  Thickness  is  Plotted  With  Different  Gear 
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CHAPTER  V 


CONCLUDING  REMARKS 


A  computer  solution  for  the  dynamic  load  in  a  pair  of  spiral  bevel 


gear  sets  was  developed  by  solving  the  equations  of  motion  for  the  pinion 


and  gear  shaft.  An  existing  finite  element  code  was  used  to  calculate 


the  combined  stiffness  of  the  contacting  pinion  and  gear  teeth  as  a 


function  of  contact  position  in  the  zone  of  action.  In  addition  to 


the  dynamic  load  analysis,  a  computer  solution  was  developed  to 


predict  the  bulk  surface  temperature,  the  flash  temperature,  and  the  film 


thickness  along  the  contact  path.  An  existing  finite  element  heat  code 


was  used  to  calculate  the  temperature  influence  coefficients  from 


which  the  bulk  surface  temperature  is  calculated.  Both  the  lubricant 


film  thickness  and  the  sliding  traction  are  calculated  from  the  recent 


findings  in  EHL  theories. 


Results  were  obtained  for  a  set  of  experimental  spiral  bevel  gears 


currently  being  tested  at  NASA  Lewis  Research  Center.  The  results  of 


dynamic  load  show  that  there  exist  numerous  peaks  in  the  variation  of 


dynamic  load  against  the  gear  shaft  speed.  These  fluctuations  correspond 


reasonably  well  with  the  critical  frequencies  of  the  system.  The  en¬ 


velope  of  the  peaks  suggests  that  the  highest  dynamic  load  occurs  some¬ 


where  near  the  critical  frequency  corresponding  to  the  rotational  mode 


oscillations  of  the  two  gears. 


Results  of  the  film  thickness  show  that  its  variation  along  the 
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contact  path  is  not  large,  and  is  caused  mainly  by  the  increase  in  the 
entrainment  velocity  and  the  change  in  the  bulk  surface  temperature.  The 
total  flash  temperature  variation  is  controlled  by  the  sliding  velocity, 
having  its  maximum  near  the  end  of  the  contact  path  where  the  transition 
from  double  to  single  teeth  pair  mesh  occurs.  Effects  of  operating 
variables  on  the  minimum  film  thickness  and  maximum  surface  temperatures 
along  the  contact  path  can  be  obtained  readily  with  this  program. 

Results  for  the  effect  of  change  in  the  ambient  viscosity  show  trends 
consistent  with  those  anticipated  from  existing  EHL  theories. 
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TABLE  II 


GEAR  DATA  AND  LUBRICANT  DATA 


Gear  Data: 

Teeth 

Pitch  Angle 
Shaft  Angle 
Spiral  Angle 
Diametral  Pitch 


Gear 

36 

71°34' 


Pinion 

12 

18°26' 


Standard  Operating  Conditions: 
Gear  RPM 
Pinion  RPM 

Load  at  Pitch  Point  N  (lb) 
Ambient  Temperature  °C  (°F) 


5000 

15000 

11800  (2660) 
37.8  (100) 


Geometric  Dimensions: 


0.1658  m 
0.0312  m 
0.04336  m 
0.1964  m 
0.2515  m 
0.0325  m 
0.09311  m 
0.1987  m 


(see  Fig.  C-l) 
(6.527  in) 
(1.230  in) 
(1.707  in) 
(7.733  in) 
(9.901  in) 
(1.280  in) 
(3.6656  in) 
(7.824  in) 


Gear  Material  Data: 
Steel: 


-» V.  V- Q 

Density  g/cin  (lb/in  ) 


7.81  (.282) 


Thermal  Conductivity  at  311°K  (100°F) 
w/m°K  (BTU/ (gee) (in) (°F)) : 

Young's  Modulus  GPfl  (psi) 

Poisson  Ratio 
Surface  Convectivity 

w/m2  °K  (BTU/ (sec) (in  ) (°F)) : 

Oil  Jet 
Oil/Air  Mist 
Air 


46.7  (.000625) 
207.  (30000000) 

0.3 


397.  (0.000135) 

19.8  (0.00000675) 
3.97(0.00000135) 


Lubricant  Data: 

Super-Re fined,  Napththenic,  Mineral-Oil 
Dynamic  Viscosity  at  311°K  (100°F) 
cp  (lb* sec/in2) 

Density  at  311°K  (100°F)  g/cm3  (lb/in3) 
Thermal  Conductivity  at  311°K  (100°F) 
w/m°K  (BTU/(sec) (in) (°F)) 
Viscosity-Pressure  •  Temperature  Relation 
u  -  U  exp(ap  +  e(^  -  =^)) 


64.7  (0.0000 0938) 
0.61(0.022) 


0.125(0.00000168) 


Pressure-Viscosity  Coefficient 
a  m2/MN  (in2/lb) 

Temperature-Viscosity  Coefficient 
S  °K  (°R) 


0.023(0.00016) 


3890  (7000) 
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NOMENCLATURE 
(For  Chapter  II) 

outer  cone  distance 

vector  from  machine  center  to  the  point  on  the  surface 
addendum 

linear  acceleration  vector 

vector  from  wheel  center  to  the  point  on  the  surface 

dedendum 

distance  from  cradle  center  to  cam  rotation  center 

nominal  wheel  diameter 

machine  coordinate 

side  dresser  offset 

blank  offset 

running  offset 

face  width 

unit  vector  along  the  gear  axis 

distance  along  tooth  axis  from  crossing  point 


inside  dresser  arm  length 
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N,  n 


n 

n. 

i 

P 

Q 

R 

R 

Y 

Yd 

S 

(S) 

T 

T 

o 

r 

u 


outside  dresser  arm  length 
velocity  ratio 

number  of  teeth  for  gear  and  pinion 
unit  normal  vector  to  the  cutter  blade 


index  interval 


unit  vector  along  the  pinion  axis 
cradle  setting  angle 

radius  from  tooth  axis 


vector  from  crossing  point  to  the  point  on  the  surface 

vector  from  the  intersection  of  plane  containing  the 
cutter  tips  and  wheel  axis  to  the  point  on  the  surface 

side  dresser  radial 


vector  from  machine  center  to  the  intersection  of  plane 
containing  the  cutter  tips  and  the  wheel  axis 

surface  expression  in  matrix  form 

cam  setting 

standard  cam  setting 

unit  vector  along  the  direction  of  dressing  diamond 
movement 

velocity  of  contact  point 


3 


X,  Y,  X 


*b 


di 


do 


gy 


py 


r,  y 


VY0 


r  ,  y 

Y  Y 


linear  velocity  vector 


velocity  of  the  contact  point  across  the  surface 


Cartesian  coordinate 


machine  center  to  back 


sliding  base 


inside  diamond  setting 


outside  diamond  setting 


gear  apex  withdraw 


pinion  apex  withdraw 


side  dresser  axial 


cam  guide  angle 


eccentric  angle 


pitch  angle  for  gear  and  pinion 


face  angle  for  gear  and  pinion 


root  angle  for  gear  and  pinion 


cam  rotation 


cam  pitch  radius 


\  I 


0 


wheel  cutter  rotation  angle  based  on  the  line  perpen¬ 
dicular  to  the  S. 


0 


d 

X 


P 


'P 


dresser  rotation  angle 

work  rotation 

radius 

shaft  angle 
running  shaft  angle 
dresser  block  angle 
outside  pressure  angle 

cradle  rotation 

spiral  angle 

angular  velocity  vector 


Subscripts  g  and  p  refer  to  gear  and  pinion,  respectively. 
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a 


i 

i 

) 

i 


F 


Fc 

Fr 

FMA 

*G 

i»  j  ,  k 

I 

J 


ks 

m 


M, 


y 

c 


Y 


r 


RBG 


NOMENCLATURE 
(For  Chapter  III) 

acceleration  vector 

external  force  vector 

tooth  contact  force 

bearing  reaction  force 

scalar  quantity  of  contact  force 

angular  momentum  about  gravity  center 

unit  vector 

central  mass  moment  of  inertia  about  the  principal  axis 
of  inertia 

equivalent  mass 

combined  tooth  contact  stiffness 

mass 

vector  of  external  moment  about  gravity  center 

position  vector  from  gravity  center  to  the  tooth 
contact  point 

position  vector  from  gravity  center  to  the  bearing 
reaction  point 

radius  of  effective  base  circle  of  gear 
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i-  / 

RBP 

radius  of  effective  base  circle  of  pinion 

T 

static  force 

> 

c 

X,  y,  Z 

translational  displacement 

.  V 

9 

rotational  displacement 

<Y*| 

! .  r 

Q 

relative  rotational  velocity  of  moving  coordinate  to 
fixed  coordinate 

0) 

nominal  shaft  velocity 

{  } 

column  matrix 

life 

1 1 

matrix  form 

V|‘‘ 

[C] 

damping  matrix 

. 

•!*:* 

[DC] 

transformation  matrix 

•  «*/ 

v>7 

[DG] 

transformation  matrix 

V  ,V 

[DK] 

transformation  matrix 

{G  <> 
ni 

normal  displacement  at  point  i  due  to  rigid  body  motion 

(K> 

stiffness  matrix 

[m] 

mass  matrix 

(S»ci} 

normal  displacement  at  point  i  due  to  tooth  contact 

■v 

force 

.  ,  •  -V  V 

-» V  V-'d  '  '  ’  ~  WTvF  V  V  V  *1’  —  7  ’ “  m 

L: ■  •  .  ■  ■  • ‘-m 
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[eu] 


transformation  matrix 


[Qk] 


transformation  matrix 


Subscripts  p  ,  g  refer  to  pinion  and  gear. 


r 
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i 

NOMENCLATURE 

| 

(For  Chapter  TV) 

a 

thermal  diffusivity 

1  % 

AMAX 

length  of  semi-major  axis  of  contact  ellipse 

I 

BMAX 

length  of  semi-minor  axis  of  contact  ellipse 

i 

E' 

effective  elastic  modulus 

W  *_•* 

G 

material  parameter 

H  *  « 

1  < ' 

G 

00 

limiting  elastic  shear  modulus  •* 

1 

h 

film  thickness 

H  <  ' 

hj 

surface  heat  transfer  coefficient  for  oil  jet  cooling 

h 

s 

surface  heat  transfer  coefficient  for  free  convective  f 

cooling  | 

ht 

surface  heat  transfer  coefficient  for  air/oil  mist 
cooling 

H  . 
mm 

dimensionless  film  thickness  of  point  contact  solution 

H  ^  0 
min,£ 

dimensionless  film  thickness  of  line  contact  solution 

1  ':V 
■  '.vl 

H 

k 

solid  material  thermal  conductivity 

9 ' 

k 

o 

oil  thermal  conductivity 

K 

ellipticity  parameter 

If  4 

l 

length  of  contact  stripe 

L 

i  ’ 

*  F  tT 

S  'V'.;:  W  V*  1 - J  ■  * %  TTf  V “  >  !  j  v  ,  ii  i  1  .  | 
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ave 


mid 
u,  V 

U 


«r 


m 


pressure 

heat  flux 

heat  flux  at  tooth  contact  point 
nominal  length,  effective  radius 
temperature 

average  temperature  over  whole  contact  ellipse 

bulk  temperature 

flash  temperature 

temperature  at  shear  plane 

speed 

speed  parameter 
sliding  velocity 
load  per  unit  length 
load  parameter 

distance  from  solid  surface  to  shear  plane 

heat  partition  coefficient,  pressure  dependent  viscosity 
parameter 

shear  rate 


dimensionless  shear  rate  (y  y  /x  ) 

O  L 


shear  rate  of  elastic  part 


shear  rate  of  viscous  part 


atmosphere  viscosity  of  lubricant 


low  shear  stress  viscosity 


shear  stress 


dimensionless  shear  stress  (t/tt) 

Li 

limiting  shear  stress 
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APPENDIX  A 

SAMPLE  DATA  FOR  MACHINE  SETTINGS 

Gear  Pinion 


Concave 

Convex 

Number  of  Teeth 

36 

12 

Face  Width 

1.00 

Shaft  Angle 

1.5708 

Outer  Cone  Distance 

3.691 

Addendum 

0.093 

0.231 

Dedendum 

0.267 

0.129 

Pitch  Angle 

1.249 

0.3217 

Face  Angle 

1.276 

0.3895 

Root  Angle 

1.1813 

0.2944 

Spiral  Angle 

0.61086 

Nominal  Wheel  Diameter 

6.0 

Machine  Center  to  Back 

0.0 

-0.067 

0.060 

Sliding  Base 

-0.016 

-0.009 

-0.046 

Blank  Offset 

0.000 

0.041 

-0.045 

Cam  Setting 

6.916 

7.089 

7.279 

Eccentric  Angle 

0.3598 

0.3616 

0.3604 

Cradle  Angle 

1.2116 

5.3354 

5.4990 

Cam  Guide  Angle 

0.0 

0.0 

0.0 

Standard  Cam  Setting 

7.00 

7.00 

7.00 

Cam  Pitch  Radius 

6.2273 

5.9567 

5.9567 

Index  Interval 

11 

11 

11 
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APPENDIX  A  (continued) 


Gear  Pinion 


Concave 

Convex 

Dresser  Block  Angle 

0.7854 

0.7854 

Outside  Pressure  Angle 

0.3490 

0.3490 

0.3482 

Side  Dresser  Radial 

3.113 

3.028 

2.959 

Outside  Dresser  Arm  Length 

1.0 

1.0 

1.0 

Outside  Diamond  Setting 

0.0 

0.0 

0.171 

Inside  Dresser  Arm  Length 

1.0 

1.0 

1.0 

Inside  Diamond  Setting 

0.223 

0.170 

0.0 

Side  Dresser  Axial 

0.0 

0.023 

0.029 

Side  Dresser  Offset 

1.00 

0.800 

1.200 

i 


1 
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APPENDIX  B 


From  the  geometry  of  Fig.  2-10,  the  motion  formulas  for  a  stand¬ 
ard  cam  rotation  e  are  obtained  as  follows: 

o 


sin  ip  +  sin  a  +  sin(e 
O  DC  O 


V 


(B-2) 

(B-3) 


(B-7) 
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where  N  =  n  cos  U >  +  An  cos(e  -  1  ) 

u  o  u  o  o 

From  (B-4)  and  (B-5) 


de 

m^2m 

dp. 

1  4.  _ £ 

dp 

0 

de 

de 

1  +  de 

de 

o  o  o  o 

sec  i n 

=  1  +  N  - - - - — 

DC  +  n  DC  +  n 
u  u 


=  D 


sec  p, 

_ i 

DC  +  n 


(B-8) 


u 


where  D  =  DC  cos  p  +  N  and  also  4^""  =  —  =  4^- 
2  de  n  de 


, ,  , ,  de  de 

dt|>  _  diji  _ o  _  dij/  _ o  __n_ 

dX  de  dX  de  de  n, 
o  o  i 


_n_  N 
n^,  D 


(B-9) 


The  term  a  = 
c 


3X2J 


used  in  determining  the  motion  parameters  is 


derived  as  follows: 


d  p  _  _n _ 1_ 

dX2  ni  D2 


D  —  -  N  — 
dX  dX 


2  , 

n 

1 

n . 

2 

ij 

D 

de  de 

) 
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•  * 

n 

2  (DC  +  n  ) 
u 

cos  \p 

fn  dN_ 

dD 

M 

n , 

3 

D  d r ' 

N  “ 

de 

iJ 

D 

o 

o 

j 

(B— 1 0 ) 


and 


dN  <  ,  <1*o  „  < 

- —  =  -  n  sin  ip  - - An  sin 

de  u  o  de  u 

o  o 


e  -  ip 
o  o 


1  - 


d* 

de 


oj 


AM 


DC  +  n 


(B-ll) 


where 


AM  =  n  sin  ip  +  An  DC  sin(e 
u  o  o 


V 


dD  _  ^2  AM 

deo  "  S  n  *2  de  ~  DC  +  n 

o  u 


(B-12) 


5E-TV  lDC  N  c“  *2  +  "» 

U 


(B-12) 


,2,  2  cos  ^9 

S-l.  „  (JL)^  -  2  {N(DC  N  tan  ^  +  AM)  -  D  AM) 

'  dX2  "i  D3  2 


(B-13) 
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APPENDIX  C 


The  details  of  the  mass,  viscous  and  stiffness  matrices  of  the 


equation  of  motion  are  shown  here. 


Number  of  contact  points  at  each  time 


AXYZCG  (I,K) 


Normal  vector  of  gear  surface  at  contact  point 


(inward  to  gear  body) 


AXYZCP  (I,K) 


Normal  vector  of  pinion  surface  at  contact 
point  (outward  to  pinion  body) 


AI  (J,K) 


Stiffness  of  contact  point  I  due  to  a  unit 


load  applied  at  contact  point  K 


XGG,  YGG,  ZGG 


Translational  displacements  of  gear  gravity 


center 


AXG,  AYG,  AZG 


Rotational  displacements  of  gear  body 


SGP ,  YGP,  ZGP 


Translational  displacements  of  pinion  gravity 


center 


AXP,  AYP ,  AZP 


Rotational  displacements  of  pinion  body 


SYZrG  (I,K) 


Position  vector  from  gear  gravity  center  to 


contact  point  K 


XYZCP  (I,K) 


Position  vector  from  pinion  gravity  center 


to  contact  point  K 


where  1=1,  2,  3  indicates  the  component  in  three-dimensional  axis. 


J  and  K  equal  one  to  MAXC.  Array  ABG  (J,I),  J  =  1  to  12,  I  -  1  to  MAXC 


was  introduced  for  convenience  in  developing  equations: 


ABG  (1,1) 


=  E  -  AI(I,K)*AXYZCG(1,K) 
K 
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ABG  (2,1) 
ABG  (3,1) 
ABG  (4,1) 

ABG  (5,1) 

ABG(6,K) 

ABG  (7,1) 
ABG  (8,1) 
ABG  (9,1) 
ABG  (10,1) 

ABG  (11,1) 

ABG  (12,1) 


*  Z  -  AI(I,K)*AXYZCG(2,K) 

K 

=  X  -  AI ( I , K) *AXYZCG (3 , K) 

K 

*  Z  -  AI(I,K) [XYZCG (2 ,K)*AXYZCG (3 ,K) 
K 

-  XYZCG(3 ,K)*AXYZCG (2 ,K) ] 


-  Z  -  AI(I,K) [AYZCG(3 ,K)*AXYZCG(1 ,K) 
K 

-  XYZCG (1 , K) *AXYZCG (3 , K) ] 


-  Z  -  AI (I,K) [XYZCG (1 ,K)*AXYZCG (2 ,K) 
K 

-  XYZCG (2 ,K)*AXYZCG (1 ,K) ] 


-  Z  AI ( I , K) *AXYZCP (1 , K) 

K 

*  Z  AI(I,K)*AXYZCP(2,K) 

K 

-  Z  AI(I,K)*AXYZCP(3,K) 

K 

-  Z  AI(I,K)[XYZCP(2,K)*AXYZCP(3,K) 
K 

-  XYZCP(3,K)*AXYZCP(2,K) ] 


-  Z  AI(I  ,K)  [XYZCP (3  ,K)*AXYZCP (1  ,K) 
K 

-  XYZCP (1 , K) *AXYZCP (3 ,  K)  ] 


=*  Z  AI  (I ,  K)  [XYZCP  (1  ,K)*AXYZCP  (2  ,K) 
K 

-  XYZCP (2 , K)*AXYZCP (1 ,K) ] 


The  summation  is  from  K  *  1 
The  force  situation 


to  K  '  KMAXC ,  and  I  =  1  to  KMAXC. 
is  shown  in  Fig.  3-1  and  Fig.  C-l. 
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For  the  gear,  the  components  of  the  contact  force  are 


AXYZCG  (1,1)  { ABG (J ,  I)  } 


+AXYZCG  (2,1)  {ABG(J, I) } 


P' 


+AXYZCG  (3,1)  (ABG(J,I)} 


J  =  1  to  12 


where  { 


denotes  column  matrix. 


D 

P 


(XGG,  YGG,  ZGG,  AXG,  AYG,  AZG , 


XGP ,  YGP,  ZGP,  AXP ,  AYP,  AZP)  and  I  from  one  to  MAXC.  The  damping  force 
at  contact  point  I  is 

the  force  terms  due  to  bearing  reaction  are 
-2KXP*XGP  +  KXP (RIP-R2P)*AZP  I 

P 

-KYP*YGP  J 

JP 

-2KZP*ZGP  +  KZP (R2P-R1P) *AXP  k 

P 


and  the  torque  due  to  these  bearing  forces  are 
- (R1P-YGP, )KZP*ZGP- (RIP- YGP0 ) 2KZP*AXP 
+(R2P+YGP)KZP-(R2P+YGP0)KZP  +  AXP  i 
+(R1P-YGP. ) KXP*XGP- (RIP-YGPo ) 2KXP*AZP 
- (R2P+YGP0 )KXP*XGP- (R2P+YGP0 ) 2KXP*AZP  k 


"^D-A184  772 
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where  YGP„  is  the  rigid  body  displacement  of  the  pinion  in  the  pinion 
fixed  coordinate. 

The  damping  force  due  to  the  bearing  is 

C’  *XGP  +  C'  *AZP  I 
xp  azp  p 

+C'  *YGP  7 

yp  P 

•  ■ 

+C'  *ZGP  +  C'  *AXP  k 
zp  axp  p 

and  the  damping  torque  is 

C"  *ZGP  +  C"  *AXP  I 


zp 


axp 


+C"  *XGP  +  C"  *AZP  k 
xp  azp  p 


Now  the  coefficients  of  damping  will  be  determined.  Let  I  *=  1 
to  MAXC  and  introduce  array  CBG(j,I),  J  =  1  to  12  for  convenience,  and 
define 


CBG  (1,1) 
CBG  (2,1) 
CBG  (3,1) 
CBG  (4,1) 
CBG  (5,1) 
CBG  (6,1) 
CBG  (7,1) 
CBG  (8,1) 
CBG  (9,1) 
CBG(10, I) 


-AXYZCG(1,I) 

-AXYZCG(2, I) 

-AXYZCG(3,I) 

- [XYZCG(2 , I)*AXYZCG (3 , I)-XYZCG (3 , I)*AXYZCG (2,1)] 
- [XYZCG ( 3 , I) *AXYZCG ( 1 , I ) -XYZCG ( 1 , I ) *AXYZCG (3 , I) ] 
- [ XYZCG (1,1) *AXYZCG (2,1) -XYZCG (2,1) *AXYZCG (1,1)] 
AXYZCP(l.I) 

AXYZCP(2,I) 

AXYZCP(3,I) 

XYZCP(2, I)*AXYZCP(3 , I)-XYZCP (3 , I)*AXYZCP(2 , I) 
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CBG  (11,1)  -  XYZCP (3,1) *AXYZCP (1,1) -XYZCP (1,1) *AXYZCP (3,1) 
CBG  (12,1)  =  XYZCP (1 , I)*AXYZCP (2 , I)-XYZCP(2 , I)*AXYZCP (1 , I) 

For  gear,  J  =*  1  to  12 

C  =  AXYZCG(l.I)  {CGB(J,I)}T  C 

-  C 

C  -  AXYZCG(2,I)  {CBG(J, I) }T  C 
cyg  c 

C  =  AXYZCG(3,I)  (CBG (j ,I)}T  C 

tig  c 

C'  "  tXYZCG(2,I)*AXYZCG(3,I)-XYZCG(3,I)*AXYZCG(2,I)] 

cxg 

(CBG(J,I)}T  Cc 

C^yg  =■  [XYZCG(3,I)*AXYZCG(1,I)-XYZCG(1,I)*AXYZCG(3,I)] 
(CBG(J,I)}T  Cc 

Cczg  *  [XYZCG(1>I)*AXYZCG(2,I)-XYZCG(2,I)*AXYZCG(1,I)] 
(CBG(J,I)}T  Cc 

C'  =  -COG 
Xg 

C’  -  -C1G  -  C2G 

yg 

C'  -  C1G* (R1G-XGG o ) ~C 2G* (R2G-XGG „ ) 

C«  «  -C1G-C2G 

zg 

Cayg  "  ~C1G* (R1G+XGG0 )+C2G* (R2G-XGG0 ) 

Czg  "  "^R1G+XGGo)ClG+(R2G-XGG0)C2G 

C"  -  -(R1G+XGGJ2  C1G-  (R2G-XGG-  ) 2  C2G 
azg 
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C"  -  (R1G+XGG  o ) C1G- (R2G-XGG  0 ) C2G 

yg 

C"  -  - (R1G+XGG 0 ) 2  C1G- (R2G-XGG  0 ) 2C2G 
ayg 

For  pinion,  J  *  1  to  12 

C  =  -AXYZCP(1,I){CBG(J,I)}T  C 
exp  c 

Ccyp  =  -AXYZCP(2,I){CBG(J,I)}T  cc 

C  *  -AXYZCP(3,I){CBG(J,I)}T  C 
cap  c 


The  torque  due  to  this  contact  force  is 

[XYZCG(2,I)*AXYZCG(3,I)-XYZCG(3, I)*AXYZCG(2,I)]  {ABG(J,I)  T 


i 


+  [XYZCG (2 , I)*AXYZCG (3 , I)-XYZXG (3 , I)*AXYZCG(2 , I) ]  {ABG(J,I)}T 


+  [XYZCG(1,I)*AXYZCG(2,I)-XYZCG(2,I)*AXYZCG(1,I) ]  {ABG(J,I)}T 


and  the  damping  torque  is 
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For  pinion  case,  the  components  of  contact  force  are 


-AXYZCP(1,I;  { ABG(J, I) } 


-AXYZCP(2,I)  (ABG(J.I)} 


D 

l  P 


T  vj, 


“s' 


-AXYZCP(3,I)  (ABG(J,I)}  (d 


I  -  1  to  MAXC  and 
J  -  1  to  12 


the  damping  force  at  contact  point  I  is 


C  pa 
exp  J  g 


P  '  P 

the  torque  due  to  this  contact  force  is 


t  i  +C 

i  ]  +c 

1  P  cyp 

ID  1  Jp  czp 
I.SI 

'  8 

) 

Id  J 

- [XYZCP (2,1) *AXYZCP (3 , I) -XYZCP (3 , I) *AXYZCP (2 , I)  (ABG(J,I)} 


vlp 


- [XYZCP (3 , I) *AXYZCP (1 , I)-XYZCP (1 , I)*AXYZCP (3 , I) ] { ABG ( J , I) } 

D  )  jp 
g  1  F 


-[XYZCP(1, I)*AXYZCP(2 , I)-XYZCP(2 ,I)*AXYZCP(1 , I) ] {ABG(J , I) } 


g 


Tt 


1  to  12  and  the  damping  torque  is 


+  C' 
cyp 

C '  -  " lXYZCP  < 2 » I) *AXYZCP (3 , 1)-XYZCP (3 , I) *AXYZCP (2 , 1) ] 

{CBG(J,I)}T  cc 

Ccyp“  " f XYZCP ( 3 ’ 1 ^  *AXYZCP  d » 1 ) -XYZCP (1 , 1) *AXYZCP (3 , 1) ] 
{CBG(J,I)}T  C 

c 

C'  -  - [XYZCP (1 , I) *AXYZCP (2 , I) -XYZCP (2 , I) *AXYZCP (1,1)] 
cz  p 

{CBG(J,I)}X  c 

C 

C'  -  -C1P  -  C2P 
xp 

C'  -  C1P* (R1P-YGP0 )  -  C2P*(R2P+YGP„) 
slz  p 

C'  -  -COP 

yp 

C'  -  -C1P  -  C2P 
zp 

C’  -  -C1P* (RIP-YGPo )  +  C2P* (R2P+YGP, ) 

*xp 

czp  “  -(Rlp-YGPo)ClP  +  (R2P+YGP0 )C2P 

Caxp"  "  (R1P"YGP<» )  2clp-  (R2P+YGP,  )  2C2P 

Cxp  "  (RIP-YGPo ) C1P- (R2P+YGPq ) C2P 

C"  -  -(R1P-YGP0)2C1P-(R2P+YGP„)2C2P 
azp  0 


•  \ 

_ 

i  ■  \ 

D 

j  +c' 

D 

I 

D  j 

P  czp 

D  1 

o' 
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APPENDIX  D 

INTERPOLATION  OF  THE  COEFFICIENTS  OF  DEFLECTION 
AND  BULK  TEMPERATURE 


(A)  Deflection  Coefficient 

There  are  thirty  nodal  points  in  the  contacting  tooth  surface. 

The  deflection  coefficient  D^_.  is  the  deflection  (normal  direction)  at 

point  i  due  to  a  unit  load  applied  at  point  j .  The  coefficient  D  was 

Aj 

obtained  from  the  finite  element  method.  Fig.  D-l  shows  these  grid  points 
in  the  axial  plane  and  also  shows  the  contact  path.  It  is  desired  to 
determine  the  deflection  D^  at  contact  point  I  in  surface  element  a 
with  four  corner  points  jl,  j2,  j3  and  j4,  due  to  a  unit  load  at  this 
point  I.  The  concept  of  the  shape  function  used  in  the  finite  element  method 
was  borrowed  here  to  solve  the  problem.  A  mapping  diagram  is  shown  in 
Fig.  D-2,  correspondent  to  Fig.  D-l.  All  the  surface  elements  in 
this  mapping  diagram  are  of  the  same  size  as  the  square  shape.  Intro¬ 
ducing  two  new  local  coordinate  variables  r  and  s,  one  can  obtain  the 
deflection  at  I  by 


Dj  -  0.25*(l-r)*(l-s)*DJ1  JX  +  0.25*(l+r)*(l-s)*DJ2 

+  0. 25*(l+r)* (l+s)*D  ,  T,  +  0.25* (1-r )* (1+s) *D  T, 

J3,J3  J4,J4 

r,  s  are  defined  as 

X  -  0.25*(l-r)*(l-s)*X  4  0. 25*(l+r)* (l-s)*X  _ 

J  Jl  JZ 

+  0. 25*(l+r )*(l+s)*Xj3  +  0.25*(l-r)*(l+s)*Xj4 
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Yj  =  0.25*(l-r)*(l-s)*YJ1  +  0. 25*(l+r)*(l-s)*YJ2 

+  0.25*(l+r)*(l+s)*Y„  +  0.25*  (1-r )*(l+s)*YT, 

Jj  J4 

which  X  and  Y  are  the  coordinates  in  the  axial  plane. 

(B)  Bulk  Temperature 

,  Cv 

The  temperature  rise  (T^)  at  contact  point  I  due  to  a  unit  heat 
flux  at  the  same  point  can  be  calculated  the  same  way  as  that  for  the 
deflection  coefficient.  If  T^  is  the  temperature  rise  at  nodal  point  i 
due  to  a  unit  heat  flux  at  nodal  point  j ,  then 

Tj  -  0.25*(1-  )*(l-s)*TJ1  J;L  +  Q.25*(l+r)*(l-s)*TJ2  J2 

+  0.25* (1+r ) * (l+s ) T J3  J3  +  0.25*(l-r)*(l-s)*Tj4  J4 

Q 

Now  the  temperature  rise  (T  _)  at  contact  point  J  due  to  a  unit 

J  ,  I 

heat  flux  at  contact  point  I  (I^J)  has  to  be  determined.  If  point  J  is 
far  away  from  point  I,  then  the  temperature  rise  at  point  J  will  not  be 
influenced  by  the  local  temperature  jump  at  point  I.  Fig.  D-3  shows 
this  situation  where  contact  I  is  in  surface  element  0  with  nodal  points 
II,  12,  13,  14  and  contact  point  J  is  in  surface  element  a  with  nodal 
point  Jl,  J2,  J3  and  J4.  The  temperature  rise  at  nodal  points  II,  12, 

13,  and  14  due  to  heat  flux  at  contact  point  J  is 

TII,J  '  0-25*<1-rJ,*(1-SJ,*IXI.Jl  +  0-25*(1+rJ)*<1-sJ>*TII>J2 

+  °-25*<l«J)*(l+sJ>*TIIjJ3  +  °-25*<l-rJ>*<l+sJ>*TII,J4 

where  II  =  II,  12,  13  and  14,  respectively,  and  r^,  s^  are  local 
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coordinates  in  surface  element  a.  Then  the  temperature  rise  at  contact 
point  I  will  be 

T5,j  -  0.25*(l-rI)»a-.t)*TI1>J  + 

+  0. 25*(l+r  )* (1+s  )*T  +  0.25*(l-r  )*(l+s  )*T 

I  II3,J  I  I  14, J 

where  and  S^.  are  local  coordinates  in  surface  element  8. 

If  contact  point  I  is  close  to  contact  point  J,  the  local 
temperature  effect  should  be  counted.  One  can  make  a  square  with  the 
same  size  as  any  other  squares  so  that  the  contact  point  I  will  be 
within  this  square  and  contact  point  J  will  become  one  of  the  four  corner 
points  of  this  square  (see  Fig.  D-4).  Then  the  temperature  at  each 
corner  point  due  to  heat  flux  at  point  J  can  be  determined  as  above. 

Let  r^.  and  s^  be  the  local  coordinates  in  the  new  surface  element,  then 

j  =  0.25*(l-rI)*(l-sI)*TJ  +  0.25*(l+rI)*(l-s];)*T1 
+  0.25*(l+rI)*(l+sI)*T2  +  0.25*(l-rI)*(l+sI)*T3 

the  order  of  T^,  T^,  and  T^  will  be  changed,  depending  on  how  this 
new  element  is  formed.  In  order  to  cover  all  the  possible  cases  where 
this  new  element  may  locate,  a  new  element  layer  was  created  to  surround 
the  original  mapping  diagram  (as  shown  in  Fig.  D-4).  The  temperature  at 
these  new  added  grid  points  due  to  heat  flux  at  any  possible  grid  point 
was  given  artificially. 


r  a/d 

jO-%1 
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